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Abstract 



We provide a systematic effective lagrangian description of the phenomenology of the 
hghtest top-partners in composite Higgs models. Our construction is based on symmetry, 
on selection rules and on plausible dynamical assumptions. The structure of the resulting 
simplified models depends on the quantum numbers of the lightest top partner and of the 
operators involved in the generation of the top Yukawa. In all cases the phenomenology is 
conveniently described by a small number of parameters, and the results of experimental 
searches are readily interpreted as a test of naturalness. We recast presently available 
experimental bounds on heavy fermions into bounds on top partners: LHC has already 
stepped well inside the natural region of parameter space. 



1 Introduction 



The exploration of the weak scale at the Large Hadron Collider is set to unveil the dynamics of 
electroweak symmetry breaking. A giant step in that direction was achieved this year with the dis- 
covery of a bosonic resonance, whose features are remarkably compatible with those of the Standard 
Model Higgs boson. Whether we like it or not, the main question now facing us concerns the role 
of naturalness in the dynamics of the newly discovered boson. Theoretically we can think of two 
broad scenarios that concretely realize naturalness: supersymmetry and compositeness. In the case 
of supersymmetry, the implications and the search strategies have been worked out in much greater 
detail than in the case of compositeness. That is explained partly by the undisputable theoretical 
appeal of supersymmetry (gauge coupling unification, connection with string theory, etc.) and partly 
by the comfort of dealing with a perturbative set up. The difficulty in dealing with strong dynamics 
has instead, and for a long time, slowed down progress in the exploration of compositeness, and, in 
particular, progress on its objective phenomenological difficulties (mostly flavor, but also precision 
tests). Interesting ideas were indeed put forward early on [1, 2], but the absence of a weakly cou- 
pled approach prevented more concrete scenarios to appear. However, in the last decade, thanks in 
particular to the holographic perspective on compositeness, semi-perturbative scenarios have been 
depicted and studied [3] ^. Even though a very compelling single model did not cross our horizon, 
we believe we have learned how to broadly depict interesting scenarios, while remaining sufficiently 
agnostic on the details (see for instance [4]). The first aspect of an interesting set up is that the 
Higgs is a pseudo-NG-boson associated with the spontaneous breakdown of an approximate global 
symmetry. The second aspect is that flavor arises from partial compositeness: the quarks and lep- 
tons acquire a mass by mixing with composite fermions. Partial compositeness, although much more 
convincing than the alternatives, does not, by itself, lead to a fully realistic flavor scenario. This is 
because of constraints from ex [5], electric dipole moments and lepton flavor violation (see Ref. [6] 
for a recent appraisal). In a realistic scenario partial compositeness should likely be supplemented 
by additional symmetries. In any case, and regardless of details, a robust feature is that the Higgs 
potential is largely determined by the dynamics associated with the top quark and the composite 
states it mixes to, the so-called top partners. That is in a sense obvious and expected, as the top 
quark, because of its large coupling, color multiplicity and numerics already contributes the leading 
quadratically divergent correction to the Higgs mass within the SM. It is nonetheless useful to have 
depicted a scenario that concretely realizes that expectation. The naturalness of electroweak sym- 
metry breaking depends then on the mass of the fermionic top-partners. That is in close analogy 
with the supersymmetric case, where naturalness is largely controlled by the mass of the bosonic 
top partners, the stops. 

The case of light stops in supersymmetry is being actively considered both theoretically and 
experimentally. One main goal is to effectively cover all regions of parameter space, without being 
swamped by the less relevant parameters. Simplified models or motivated assumptions like "natural 
susy" [7] offer a convenient way to achieve that goal. In the reduced parameter space (featuring 
stop mass parameters and possibly the gluino or lightest neutralino mass), the constraints from of 
experimental searches offer a direct and largely model independent appraisal of naturalness. The 

^By semi-perturbative here we mean that in these models, typified by warped compactifications, there exists a 
sufficiently interesting subset of questions, even involving physics at energies well above the weak scale, that can be 
addressed using perturbation theory. 
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Figure 1: Schematic picture of the spectrum. 

goal of this paper is to provide a similar simplified approach to describe the results of experimental 
searches for top partners. We will focus on the composite Higgs scenario based on the minimal coset 
SO{5)/ S0{4:). The basic simplifying assumption is that the spectrum has the structure depicted 
in figure 1, where one SO (4) multiplet of colored Dirac fermions ^ is parametrically lighter than 
the other states. As already illustrated in Ref. [8] for the case of bosonic resonances, in that limit 
one expects the dynamics of ^' to be described by a weakly coupled effective lagrangian. Therefore 
the simplified model, at leading order in an expansion in loops and derivatives, can be consistently 
described by a finite number of parameters. Moreover symmetry and selection rules, via the Callan- 
Coleman-Wess-Zumino (CCWZ) [9] construction, reduce the number of relevant parameters. It is 
obviously understood that the limiting situation presented by the simplified model is not expected to 
be precisely realized in a realistic scenario. However, a realistic situation where the splitting with the 
next-to-lightest multiplet is of the order is qualitatively already well described by the simplified 
model. Only if the splitting were parametrically smaller than m-q, would there be dramatic changes. 
We should also stress that our models are truly minimal, in that they do not even possess sufficient 
structure (states and couplings) to make the Higgs potential calculable. In principle we could add 
that structure. For instance by uplifting our multiplet ^' to a full split 50(5) multiplet, like in a 
two site model, we could make the Higgs potential only logarithmically divergent, thus controlling 
its size in leading log approximation, and making the rough connection between m>ii and naturalness 
more explicit along the lines of [10]. We could even go as far as making the one loop Higgs potential 
finite with a three site model [11, 12], or by imposing phenomenological Weinberg sum-rules [13]. 
However in these less minimal models the first signals at the LHC would still be dominated by the 
lightest SO{4:) multiplet, whatever it may be. The point is that while the contribution of the heavier 
multiplets does not decouple when focussing on a UV sensitive quantity like the Higgs potential, it 
does decouple when considering the near threshold production of the lightest states. For the purpose 
of presenting the results of the LHC searches in an eloquent way, the simplified model is clearly the 
way to go. There already exists a literature on simplified top partner models in generic composite 
Higgs scenarios [14, 15, 16], where the role of symmetry is not fully exploited. Focussing on the 
minimal composite Higgs model based on SO{5)/ SO{4:), our paper aims at developing a systematic 
approach where all possible top partner models are constructed purely on the basis of symmetry 
and selection rules. 

In the end we shall derive exclusion plots in a reduced parameter space, which in general involves 
the mass and couplings of the top-partner ^. Now, even though these are not the parameters of a 
fundamental model, given their overall size, we can roughly estimate how natural the Higgs sector 
is expected to be. We can then read the results of searches as a test of the notion of naturalness. To 
make that connection, even if qualitative, we must specifiy the dynamics that gives rise to the top 
Yukawa. As discussed in [17], there are several options, each leading to a different structure of the 
Higgs potential and thus to a different level of tuning. The common feature of all scenarios is that 



2 



the top partners need to be light for a reasonably natural theory, the way the tuning scales with the 
top-partners' mass is instead different in each case. In this paper we focus on the possibility that 
the right handed top quark tR is a SO{A) singlet belonging to the strong sector, therefore the top 
Yukawa simply arises from an 50(5) breaking perturbation of the form 

XLQLOR + h.C. (1.1) 

Here Or is a composite operator, which in the low energy theory maps to Htji, thus giving rise to 
a top Yukawa coupling yt ~ Xl- The operator Or however also interpolates in general for massive 
states, the top partners. Now, from simple power counting, and also from explicit constructions [11], 
at leading order in the breaking parameter Xl we expect the Higgs potential to have the form 

+2/2 + 3!/^ 



VW = -^\-'^' + ^i, + T/-h + ---}- (1-2) 



where a,b,c, . . . are coefficients expected to be 0(1), / is the decay constant of the cr-model, while 
m* broadly indicates the mass scale of the top partners. Then, since is, ideally, the lightest top- 
partner we have niys, ^ m*. Given m* and /, the measured values v = (h) = 246 GeV and m/j = 125 
GeV, may require a tuning of a and b below their expected 0(1) size. More explicitly one finds 

,= -i4vr2^f430GeVy ^^3^ 

and, defining the top-partner coupling as = m^./ f according to Ref. [4], 

k=<^_^^±_ (1.4) 
35(2 

By these equations we deduce that in the most natural scenario the top partners should not only 
be light (say below a TeV) but also not too strongly coupled. While of course the whole discussion 
is very qualitative, we still believe eqs. (1.3)-(1.4) give a valid rule of thumb for where the top 
partners should best be found. It is with eqs. (1.3)-(1.4) in mind that one should interpret the 
results of the searches for top partners. Notice that while naturalness favors sub- TeV fermionic 
resonances, electroweak precision constraints favor instead bosonic resonances above 2-3 TeV. A 
technically natural and viable model should therefore be more complex than a generic composite 
model described by a single scale. This situation closely resembles that of supersymmetric models, 
where the light squark families and the gluinos are pushed up by direct searches, while technical 
naturalness demands the stops to be as light as possible. 

This paper is organized as follows. In Section 2 we discuss the structure of the models and their 
main features such as the mass spectrum and the couplings of the top partners. Then, in Section 3 
we turn to analyse the phenomenology of the top partners, their production mechanisms and decay 
channels, highlighting the most relevant channels to focus LHC searches on. The bounds on the 
model parameters are derived in Section 4, using the LHC data available at present Finally, our 
concluding remarks are collected in Section 5. 



^ While this work was being completed ATLAS [18] and CMS [19] presented dedicated searches for top partners, 
which we did not include in our analysis. From a preliminary investigation we expect mild changes in our results from 
these new data because both the ATLAS and the CMS searches are optimized to detect pair production. As we will 
discuss in the conclusions, a radical improvement of the bounds could perhaps be achieved, with the present energy 
and luminosity, but only with searches dedicated to single production. 
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2 The Models 



Our first goal is to develop a simplified description of the top partners, suited for studying the 
phenomenology of their production at the LHC. These simplified models should capture the robust 
features of more complete explicit constructions'^ or, better, of a putative general class of underlying 
theories. In particular, robust, and crucial, features are the pNGB nature of the Higgs and the 
selection rules associated with the small breaking of the corresponding global symmetry. We will see 
below that these features strongly constraints the structure of the spectrum and of the couplings of 
the top partners, similarly to what was found in Ref. [11] for the case of partial compositeness. 

We thus assume that the Higgs is the pNGB of the minimal coset SO(5)/SO(4) and construct 
Lagrangians that respect the non- linearly realized SO (5) invariance. We follow the standard CCWZ 
construction [9], whose detailed formulation for our coset is described in Appendix A. The central 
objects are the Goldstone boson 5x5 matrix U and the and symbols constructed out of U and 
its derivative. The top partner field ^ has definite transformation properties under the unbroken 
SO (4) group. We will consider two cases, ^ transforming in the r,j, = 4 or rij, = 1 of SO (4). 

In our construction the right-handed top quark tji emerges as a chiral bound state of the strong 
dynamics, tji must thus belong to a complete multiplet of the unbroken subgroup S0(4), and, given 
we do not want extra massless states, it must be a singlet. That does not yet fully specify its 
quantum numbers. This is because, in order to reproduce the correct hypercharge, one must enlarge 
the global symmetry by including an extra unbroken \J{1)-^ factor and define the hypercharge as 
y = TJ, + X, where is the third SU(2)j:j generator of S0(5).'' Therefore the coset is actually 
S0(5) X U(l)j^/S0(4) X U(l)^, tR has X charge equal to 2/3 while the Higgs is X neutral (its 
hypercharge coincides with its charge). 

A second assumption concerns the coupling of the elementary fields, i.e. the SM gauge fields 
and Bfj^ and the elementary left-handed doublet ql = [tLibL)-, to the strong sector ^ . The EW bosons 
are coupled by gauging the SM subgroup of S0(5) x U(l)j^. The is assumed to be coupled linearly 
to the strong sector, following the hypothesis of partial compositeness [2]. In the UV Lagrangian 
this coupling has therefore the form 

£^1 = yrLKio<^'° + h.c. = y [Ql)i^ O'o + h.c. , (2.1) 

where O is an operator of the strong sector that transforms in some representation of SO (5) x 
U(l)jf . The choice of is, to some extent, free. Minimality, and the aim of reproducing explicit 
models considered in the literature, led us to consider two cases: tq = 52/3 ^i^d = ^^2/^- ^ 
Notice that the 11(1)^^ charge of the operators must be equal to the one of the in order for the 
top mass to be generated after EWSB. In total, depending on whether the top partners will be in 
the 42/3 or in the I2/3 of the unbroken 50(4), we will discuss four models named M45, M4i4 and 
MI5, Ml 14 respectively. The classification of the various models is summarized in Table 1. 

■^See [17] for a complete calculable model with totally composite Ir, analogous holographic 5d models could be 
formulated following the approach of Ref. [3] . 

''See Appendix A for the explicit form of the generators. 

^The light quark families and the leptons will not be considered here because their couplings are most likely very 
weak. 

^ Another possible option considered in the literature is ro = 4i/6. However this option is not available once tn 
is chosen to be a SO{A) singlet: the top would not acquire a mass. It should also be remarked that, regardless of the 
nature of tn, ro ~ 4i/a is disfavored when considering dangerous tree level corrections to the Zbb vertex [20, 21]. 
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ro = ^2/3 


ro = 142/3 


'^"I- = 42/3 


M45 


M4i4 


= I2/3 


MI5 


MI14 



Table 1: The nomenclature of the four models considered in the present paper, defined by the choices of the 
representations r^,ra. 



The explict breakdown of SO{5) due to y in eq. (2.1) gives rise to a leading contribution to the 
Higgs potential V{h). However, in order to be able to tune the Higgs vacuum expectation value v to 
be much smaller that its natural scale /, one may need to tune among themselves contributions to 
V{h) with a different functional dependence on h/f. In the case of tq = 142/3, the top Yukawa seed 
y itself gives rise to two independent structures, whose coefficients can be so tuned that v/ f <C 1. 
On the other hand, in the case of tq = 52/3, the leading contribution to the potential consist of 
just one structure oc siv? h/ f cos^ h/ f , with well defined, non-tunable, minima and maxima. In 
the latter case then, in order to achieve u <C /, one should assume there exists an additional of 
SO{h) breaking coupling whose contribution to the potential competes with that of the top. If this 
additional coupling does not involve the SM fields, which seems resonable, then its contribution to 
V will arise at tree level. In order not to outcompete the top contribution, which arises at loop 
level, then this coupling should be so suppressed that its relative impact on strong sector quantities 
is of order 0(y^/167r^). The latter should be compared to the effects of relative size (y/g^si)'^ induced 
at tree level by the mixing in eq. (2.1) and accounted for in this paper. We conclude that, even 
when an extra 50(5) breaking coupling is needed, it is not likely to affect the phenomenology of top 
partners in a quantitatively significant way. 

Now back to the top partners. Our choices of their quantum numbers correspond to those 
obtained in explicit constructions. However our choice could also be motivated on general grounds 
by noticing the operators O interpolate for particles with the corresponding quantum numbers. By 
decomposing O under the unbroken SO (4) we obtain, respectively, 62/3 = 42/3 + I2/3 and 142/3 = 
42/3 + 12/3 + 92/3- In both cases we expect to find a 42/3 and/or a I2/3 in the low-energy spectrum. 
It could be also interesting to study top partners in the 92/3, but this goes beyond the scope of the 
present paper. 

The coupling of eq. (2.1) breaks the S0(5) x U(l)j^ symmetry explicitly, but it must of course 
respect the SM group. This fixes unambiguously the form of the tensor A and thus of the embeddings, 
{Ql)io — ^aioQjji °f the elementary qi in S0(5) x multiplets. For the 5 and the 14, 

respectively the fundamental and the two-indices symmetric traceless tensor, we have 
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(2.2) 



Though explicitly broken, the S0(5) x U(l)j^ group still gives strong constraints on our theory. 
Indeed the elementary-composite interactions of eq. (2.1) formally respect the symmetry provided 
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we formally assign suitable transformation properties to the embeddings. Under g G SO (5) we have 

{Ql), ^ gf {Ql)v , {Ql^)ij ^ gfaf {Ql^)vj' - (2-3) 

while the U(l)jf charge is equal to 2/3 in both cases. We will have to take into account this symmetry 
in our constructions. 

2.1 Effective Lagrangians 

Based on the symmetry principles specified above we aim at building phenomenological effective 
Lagrangians for the qi, the composite t/j and the lightest top partner states The basic idea is 
that our Lagrangians emerge from a "complete" theory by integrating out the heavier resonances 
in the strong sector. We thus need to rely on some qualitative description of the dynamics in order 
to estimate the importance of the various effective operators. We follow the "SILH" approach of 
Ref. [4] and characterize the heavy resonances in terms of a single mass scale m* and of a single 
coupling = m.^/ f. As we already suggested in the introduction, parametrizing the strong sector 
in terms of a single scale is probably insufficient: a 125 GeV Higgs suggests that the mass scale of 
the fermionic resonances should be slightly lower than that of the vectors. For our purposes the 
relevant scale m* should then be identified with the mass scale of the fermionic sector. We thus 
adopt the following power-counting rule 

where 11 = n^' "'^ denotes the canonically normalized four real Higgs field components and / is 
the Goldstone decay constant. Notice the presence of the coupling y that accompanies (due to 
eq. (2.1)) each insertion of the elementary q^. Analogously the operators involving the SM gauge 
fields, omitted for shortness from eq. (2.4), should be weighted by gsM/n^*- The tR is completely 
composite and therefore it obeys the same power-counting rule as the top partner field ^. 

Two terms in our effective Lagrangian will violate the power-counting. One is the kinetic term 
of the elementary fields, which we take to be canonical, while eq. (2.4) would assign it a smaller 
coefficient, (y/g*)'^ in the case of fermions and (g/g*)'^ in the case of gauge fields. This is because 
the elementary field kinetic term does not emerge from the strong sector, it was already present in 
the UV Lagrangian with 0(1) coefficient. Indeed it is precisely because their kinetic coefficient is 
bigger than what established in eq. (2.4), that the elementary fields have a coupling weaker than 
5*. The other term violating power-counting is the mass of the top partners, which we denote by 
Mijj. We assume Mij, < in order to justify the construction of an effective theory in which only 
the top partners are retained while the other resonances are integrated out. The ratio M^/m* is 
our expansion parameter. We will therefore obtain accurate results only in the presence of a large 
separation, <C m^, among the lightest state and the other resonances However already for 

^ An organizing principle, termed partial UV completion (PUV), to consistently construct an effective lagrangian 
for a parametrically light resonance was proposed in Ref. [8] . There the focus was on the more involved case of vector 
and scalar resonances. According to PUV, the couplings involving the lighter resonance should roughly saturate the 
strength when extrapolated at the scale m«. We refer to Ref. [8] for a more detailed discussion. The effective 
lagrangians we construct in this paper automatically satisfy PUV in the range of parameters suggested by the power 
counting rule in eq. (2.4). 
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a moderate separation, Mq, < m=K, or even extrapolating towards Mij, ~ m^, our models should 
provide a valid qualitative description of the relevant physics. Nevertheless for a more careful study 
of the case of small separation our setup should be generalized by incorporating more resonances in 
the effective theory. 



2.1.1 Top partners in the fourplet 



First we consider models M45 and M4i4, in which the top partners are in the 42/3. In this case the 
top partner field is 



and it transforms, following CCWZ, as 



1 

71 



iT + 1X2/3 
\-T + X2/3/ 



(2.5) 



J 1 



(2.6) 



under a generic element g of S0(5). The 4x4 matrix h is defined by eq.s (A. 7) and (A. 8) and 
provides a non-linear representation of the full S0(5). The four ^ components decompose into two 
SM doublets {T,B) and {X^/3, X2/3) of hypercharge 1/6 and 7/6 respectively. The first doublet has 
therefore the same quantum numbers as the (tL;^L) doublet while the second one contains a state 
of exotic charge 5/3 plus another top-like quark X2/3. 

When the ql is embedded in the 62/3, i-e. in model M45, the leading order Lagrangian is 



M45 



i qlIP Ql + i tRlp tR + i ^{Ip + i^)^ - M^^^ 



+ 



ci (*r) . 7^^;. tR + yf (QlYUii ¥^ + yc2f [QlYUi^tR + h.c 



(2.7) 



where ci^2 are coefficients expected to be of order 1. Notice the presence of the ^ = e^^y^ term which 
accompanies the derivative of the top partner field: it reconstructs the CCWZ covariant derivative 
defined in eq. (A. 24) and is essential to respect 50(5). In the second line of the equation above 
we find, first of all, a direct interaction, not mediated by the coupling y, among the composite Ir 
and the top partners. This term is entirely generated by the strong sector and would have been 
suppressed in the case of partial tR compositeness. It delivers, looking at the explicit form of in 
eq. (A. 19), couplings involving the top, the partners and the SM gauge fields. These will play an 
important role in the single production and in the decay of the top partners. The last two terms give 
rise, in particular, to the top quark mass but also to trilinear couplings contributing to the single 
production of top partners. Notice that the indices of the embedding Q\ can not be contracted 
directly with those of ^ because they live in differnt spaces. The embeddings transform linearly 
under S0(5) as reported in eq. (2.3) while ^ transforms under the non-linear representation h. For 
this reason one insertion of the Goldstone matrix, transforming according to eq. (A. 7), is needed. 

For brevity we omitted from eq. (2.7) the kinetic term of the gauge fields and of the Goldstone 
Higgs, the latter is given for reference in eq. (A. 21). Moreover we have not yet specified the covariant 
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derivatives associated with the SM gauge group, these are obviously given by 



D^,qL = {d^-igW'J^-i^-g'B^-igsG^qL, (2.8) 
D^,tR = (d^-ilg'B^-igsGAtR, (2.9) 



3^ 

D^,^ = {d^-i^^g'B^-igsG^^ . (2.10) 

where g,g' and gs are the SU(2)l x U(l)y and SU(3)c gauge couphngs. We remind the reader that 
the top partners form a cofor triplet, hence the gluon in the above equation. 

The Lagrangian is very similar for model M4i4, where the ql is embedded in the symmetric 
traceless Q]^. We have 



^M4i4 ^ iqLlJ)qL + itRlJ)tR + i-^{lJ) + i^)^ - M^,-^^ 



+ 



; ci {^r)^ rd"^ tR + yf {Qiy'UiiUj, ^^R + ^/ {QiyUi,Uj5tR + h.c. 

ll4 



(2.111 



notice that the two indices of Qjf" are symmetric and therefore the term that mixes it with ^' is 
unique. The factor ^ introduced in the last term is merely conventional. 

In both models M45 and M4i4 the leading order Lagrangian contains four parameters, {M^, 
Hi ci, C2}, on top of the Goldstone decay constant /. One parameter will however have to be fixed 
to reproduce the correct top mass, while the remaining three parameters could be traded for two 
physical masses, for instance Tnx^^^ and m^, and the coupling ci. It will often be convenient to 
associate the mass with a coupling g^ 

g. = j.. 

We will see below that ci x gq, controls the strength of the interactions between the top partners 
and the Goldstone bosons at energy ~ Mq,. In particular it controls the on-shell couplings relevant 
for single production and for two body decays. Notice that, as a function of energy, the effective 
strength of this trilinear interaction is instead ~ ciE/f. For ci = 0(1), as suggested by power 
counting, the effective coupling is of order g^ = m^/ f at the energy scale of the heavier resonances, 
in accord with the principle of partial UV completion proposed in Ref. [8]. Power counting and 
partial UV completion then equivalently imply ci = 0(1) and therefore cig^ < g^. This result 
obviously follows from the fact that the Higgs is a derivatively coupled pNGB. It would be lost if 
the Higgs was instead treated as a generic resonance. In the latter case the expected coupling would 
be independent of the mass and it would be larger, of order g^,. Moreover notice that, although on 
shell it leads to an effective Yukawa vertex, the interaction associated with ci does not affect the 
spectrum when H acquires a vacuuum expectation value. That again would not be true if we did 
not account for the pNGB nature of H. The pNGB nature of H is not accounted for in the first 
thorough work on simplified top partner models [15] and in the following studies (see in particular 
[16, 22]). 

Notice that, a priori, one of the four parameters describing the simplified model could be complex. 
This is because we have at our disposal only 3 chiral rotations to eliminate the phases from the 
Lagrangians (2.7) and (2.11). Nevertheless we are entitled to keep all the parameters real if we 
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demand the strong sector respects a CP symmetry defined in Appendix A. It is easy to check that 
CP requires the non-derivative couplings to be real while the coefficient of the term involving to 
dfj, must be purely imaginary. CP conservations is an additional hypothesis of our construction, 
however the broad phenomenology does not significantly depend on it. 



2.1.2 Top partners in the singlet 



The Lagrangian is even simpler if the top partners are in the 12/3- In this case we only have one 
exotic top-like state which we denote as T. For the two models, MI5 and MI14 that we aim to 
consider the Lagrangian reads, respectively 



MI5 



: QL ilP QL + tR ilP tR + i-^ilp-^ - M^^^f 



+ 



Mil 



yf {QIYUi^^r + yc2f (QlYUi^tR + h.c. 

QL ilP QL + tR i0 tR + i^ilp-^ - M*^^' 
"|/ {qI^Ui^Uj^^r + ^/ {Q\yUi,Uj,tR + h.c. 



(2.12) 



Notice that we could have also written a direct mixing among tR and ^ because the two fields 
now have identical quantum numbers. However this mixing can obviously be removed by a field 
redefinition. Models MI5 and MI14, apart from /, contain three parameters, {M^, y, C2}, one of 
which must again be fixed to reproduce the top mass. We are left with two free parameters that 
correspond to the coupling C2 and to the mass rrij, of the partners. Notice that in this case all the 
parameters can be made real by chiral rotations without need of imposing the CP symmetry. The 
latter symmetry is automatically respected in models MI5 and Ml 14. 

In order to complete the definition of our models let us discuss the theoretically expected size 
of their parameters. From the discussion in the introduction and from experience with concrete 
models, one can reasonably argue that the favorite range for mip is between 500 GeV and 1.5 TeV, 
while g-qi is favored in the range 1 ^ (7<i' ^ 3. It is also worth recalling the favorite range of the decay 
constant / = mij- / g-;^ , which is conveniently traded for the parameter ^ defined in Ref . [3] 



(2.13) 



where v = 2mw / g = 246 GeV is the EWSB scale. Since ^ controls the deviation from the SM 
at low energies it cannot be too large. Electroweak precision tests suggest ^ ~ 0.2 or ~ 0.1, 
which corresponds to / ~ 500 GeV or / ~ 800 GeV. Smaller values of ^ would of course require 
more tuning. Finally, the strength of the elementary-composite coupling y is fixed by the need of 
reproducing the correct mass of the top quark. We will see in the following section that this implies 
y ~ = 1. 



2.2 A first look at the models 

Now that the models are defined let us start discussing their implications. The simplest aspects will 
be examined in the present section while a more detailed analysis of their phenomenology will be 
postponed to the following one. 
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2.2.1 The Spectrum 

We start from model M45 and we first focus on the fermionic spectrum. The mass-matrix after 
EWSB is easily computed form eqs. (2.7) and (2.2) by using the explicit form of U on the Higgs 
VEV obtained from eq. (A. 12). By restricting to the sector of 2/3-charged states we find 



(2.14) 



where e = (/i) // is defined as the ratio among the VEV of the Higgs field and the Goldstone decay 
constant. The relation among (h) and the EWSB scale is reported in eq. (A. 23), from which we 
derive 

C=^ = sin2e. (2.15) 

We immediately notice a remarkable feature of the mass- matrix (2.14): only the first line, i.e. the 
terms which involve the t_L, is sensitive to EWSB while the rest of the matrix remains unperturbed. 
This is due to the fact that the Higgs is a pNGB and therefore its non-derivative interactions can 
only originate from the breaking of the Goldstone symmetry S0(5). The S0(5) invariant terms just 
produce derivative couplings of the Higgs and therefore they cannot contribute to the mass-matrix. 
Since the Goldstone symmetry is broken exclusively by the terms involving the elementary qi it is 
obvious that the mass-matrix must have the form of eq. (2.14). Notice that this structure would have 
been lost if we had not taken into account the pNGB nature of the Higgs. Indeed if we had treated 
the Higgs as a generic composite S0(4) fourplet, Yukawa-like couplings of order g^, and involving 
and ^ would have been allowed. After EWSB those terms would have given rise to (2, 1) and (3, 1) 
mass matrix entries of order g^^v. 

The peculiar structure of the mass- matrix has an interesting consequence. It implies that only 
one linear combination of T and X2/3, with coefficients proportional to the (1,2) and (1,3) entries, 
mixes with the ql, while the orthogonal combination does not mix either with the ql or with any 
other state. Explicitly, the two combinations are 

1 



r 



cos^ I + sin'* I 
1 



cos 



-T + sin -X2/3 
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2 e Y 
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sm -1 
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(2.16) 



After this field redefinition the mass-matrix becomes block-diagonal 
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(2.17) 



so that the state is already a mass eigenstate with mass "1X2/3 ~ -^"t- But the spectrum also 



contains a second particle with exactly the same mass. Indeed the X5/3 cannot mix because it is the 
only state with exotic charge and therefore it maintains the mass w-Xg/g = it had before EWSB. 
The X2/3 and the X5/3 are thus exactly degenerate. This remarkable property is due to the pNGB 
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Figure 2: The typical spectrum of tlie top partners. 



nature of the Higgs and it would be generically violated, as previously discussed, if this assumption 
was relaxed. This result also depends on tji being a composite singlet. If t/j was instead a partially 
composite state mixing to a non-trivial representation of 5*0 (5) (for instance a 5) there would be 
additional entries in the mass matrix. In a sense our result depends on y being the only relevant 
parameter that breaks SO (5) explicitly. 

Once the mass-matrix has been put in the block-diagonal form of eq. (2.17) it is straightforward 
to diagonalize it and to obtain exact formulae for the rotation matrices and for the masses of the 
top and of the T partner. However the resulting expressions are rather involved and we just report 
here approximate expressions for the masses. We have 



Prom the above equation we obtain the correct order of magnitude for the top mass if, as anticipated, 
y ~ yt and (7<i< > 1 • In this region of the parameter space the corrections to the approximate formulae 
are rather small, being suppressed by both a factor ip' /g% (which is preferentially smaller than one) 
and by ^ <C 1. However we will consider departures from this theoretically expected region and 
therefore we will need to use the exact formulae in the following sections. 

Similarly we can study the sector of —1/3 charge states. It contains a massless 6^, because we 
are not including the hpc in our model, plus the heavy B particle with a mass 



This formula is exact and shows that the bottom sector does not receive, in this model, any con- 
tribution from EWSB. By comparing the equation above with the previous one we find that the 

*The top partner's spectrum with partially composite Ir has been worked out in Ref. [11, 10]. 




(2.18) 




(2.19) 
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splitting among T and B is typically small 



2 2 ^2 5>Ii + (1 ^2)^ . 2 

2 (ffi, + r) 



(2.20) 



and positive in the preferred region gij, > y, although there are points in the parameter space where 
the ordering rriT > ms can occur. The splitting among the two doublets is instead always positive, 
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y f . The typical spectrum of the top partners that we have in our model is depicted 

in figure 2. 

The situation is not much different in model M4i4. The mass-matrix for charge 2/3 states has 
again the form of eq. (2.14) 
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(2.211 



and again it can be put in a block-diagonal form by a rotation among the T and the X2/3 similar 
to the one in eq. (2.16). Therefore also in model M4i4 the physical X2/3 has mass Mij, and it is 
degenerate with the X5/3. The approximate top and T mass are given in this case by 

mt — 
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Similarly we can compute the mass of the B partner and we find 
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(2.22) 



(2.23) 



In this case, differently from model M45 (see eq. (2.19)), the mass of the B is sensitive to EWSB. 
Apart from this little difference the spectrum is very similar to the one of model M4i4 described in 
figure 2. 

The models with the singlet are much simpler because there is only one exotic state. The mass 
matrices read: 
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for models MI5 and MI14 respectively. The mass eigenvalues for model MI5 are 

C2y/ 



mt 



rrif ~ 



V2 



sme 



1 + 

1 , • 2 ' 

1 + sm e 
4^1 



(2.24) 
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(2.26) 



12 



For model MI14 instead we have 
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(2.27) 



As one can see from the last expressions the mass of the T receives positive contributions proportional 
to and hence for a fixed mass of the T, y must be limited from above. Unlike the models with 
fourplet partners, in the singlet case y completely controls the couplings of the T with the top and 
bottom quarks (see Sec. 3.2). Therefore one can expect that for a given rrif there exists a maximal 
allowed coupling of the SM particles with the top partner and hence for small masses the single 
production of T is suppressed. In addition small values of become unnatural since they require 
very small y together with a very large C2 needed to recover correct top mass. By minimizing the 
largest eigenvalue of the mass matrix with respect to for fixed y and / one can find a minimal 
allowed mass of the T which is given by 



min, MI5 

m~ 

T 



min, MI14 

m~ 

T 



mt + -^y/sine, 
m + ^^y/sin2e, 



for the models MI5 and MI14 respectively. The bound given in eq. (2.25 
plots in the following. 



(2.28) 



will affect the exclusion 



2.2.2 Trilinear Couplings 

Other interesting qualitative aspects of our models are discovered by inspecting the explicit form 
of the Lagrangians in unitary gauge. These are reported in Appendix B, and are written in the 
"original" field basis used to define the Lagrangians in eq.s (2.5, 2.7, 2.11, 2.12), i.e. before the 
rotation to the mass eigenstates. Appendix B contains, for reference, the complete Lagrangian 
including all the non-linear and the derivative Higgs interactions. However the coupling that are 
relevant to the present discussion are the trilinears involving the gauge fields and the Higgs in the 
models M45 and M4i4, reported in eq. (B.l), (B.2), (B.3) and (B.4). 

The first remarkable feature of eq. (B.2) is that the Z boson couplings with the B is completely 
standard: it is not modified by EWSB effects and coincides with the familiar SM expression gz = 
g/cw{T^ — Q). In particular it coincides with the ZhLhi coupling, involving the elementary 6^, 
because and B have the same SU{2) x U{1) quantum numbers. The Z-boson coupling to charge 
— 1/3 quarks is therefore proportional to the identity matrix. Consequently the Z interactions remain 
diagonal and canonical even after rotating to the mass eigenbasis. In particular, in the charge —1/3 
sector, there will not be a neutral current vertex of the form B — ?■ Zb. 

This property is due to an accidental parity, Plr, defined in Ref. [8] as the exchange of the Left 
and the Right S0(4) generators. This symmetry is an element of 0(4) and it acts on the top partner 
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fourplet of eq. (2.5) and on the Higgs field 11 through the 4x4 matrix 
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(2.29) 



The action of Plr is readily uplifted to 0(5) with the 5x5 matrix = diag(— 1, —1, —1, 1, 1). 
We see that Pm is not broken by the Higgs VEV, which only appears in the last component of the 
n vector. In Ref. [20], it was shown that P^fi invariance protects the Z couplings from tree-level 
corrections at zero momentum transfer. That case applied to b quarks, but the statement generalizes 
straighforwardly: if all the particles with a given charge have the same Plr, then, at tree level in the 
weak interactions, the neutral current vertices in that charge sector are canonical and, in particular, 
diagonal. 

The Lagrangians (2.7) and (2.11) are approximately Plr invariant, with the breaking coming 
only from the weak gauge couplings and from the weak mixing y between elementary and composites. 
However at tree level, for which case the elementary fields can be treated as external spectators, even 
this weak breaking is ineffective in the charge —1/3 and 5/3 sectors. Notice indeed that according 
to eqs. (2.5,2.29), under Plr, B and are odd, while T and —X2/3 are interchanged. Then, 

inspection of the embedding in eq. (2.2) shows that while we cannot assign a consiste Plr to II, 
we can instead assign negative Plr to Bl- At tree level, tL will not affect processes involving only 
quarks with charge —1/3 and 5/3, and therefore the associated explicit breaking of Plr will be 
ineffective. Analogously the breaking in the gauge sector is seen not to matter at tree level. For a 
detailed discussion we refer the reader to section 2.4 of Ref. [21]. This explains the result previously 
mentioned for the (6, B) sector and also predicts that the coupling of the X5/3 must be canonical as 
well. This is indeed what we see in eq. (B.2). 

The same argument applies to Tr and tpt in the singlet models MI5 and MI14 . The Z-vertex 
of those states is not modified, in particular there is no inZTn vertex and the production/decay 
with Z is always controlled by left-handed coupling. On the other hand, for Tl and tL the argument 
does not apply, regardless of Plr, given Tl and II do not have the same SU{2) x U{1) quantum 
numbers. 

Another interesting property concerns the W couplings of the B with the charge 2/3 states. We 
see in eq. (B.2) that the linear combination of the T and the X2/3 that couples with the B is exactly 
orthogonal to the physical (mass-eigenstate) X2/3' field defined in eq. (2.16) for model M45. Only the 
T' couples to W B, leading after the second rotation to transitions among the physical t, T and b, B. 
Such couplings are instead absent for the physical X2/3 which therefore, cannot decay to Wb. This 
feature is not, for what we can say, the result of a symmetry, but rather an accidental feature of model 
M45. In model M4i4 instead the coupling is allowed because the physical X2/3 (see the mass-matrix 
in eq. (2.21)) is not anymore orthogonal to the combination that couples to W B. Nevertheless the 
X2/S-B coupling is suppressed by (/i)^//^ and therefore the decay X2/3 — )■ Wb, though allowed in 
principle, is phenomenologically irrelevant as we will discuss in the following section. 

A final comment concerns the couplings of the physical Higgs field p. The couplings from the 
strong sector are only due to the term in eq. (B.l) and are purely derivative. Therefore, because 
of charge conservation, they cannot involve the B partner. Higgs couplings in the —1/3 charge sector 
could only emerge from the elementary-composite mixings. However they are accidentally absent 



14 



in model M45, as shown by eq. (B.3). Therefore the decay of the B to the Higgs is absent in this 
modeh In model M4i4, on the contrary, this decay is ahowed through the vertex in eq. (B.4). 

3 Top Partners Phenomenology 

Let us now turn to discuss the main production mechanisms and decay channels of the top partners 
in the models under consideration. We will first of all, in sect. 3.1, describe how the cross-sections of 
the production processes and the partial decay widths can be conveniently parametrized analytically 
in terms of few universal functions, extracted from the Monte Carlo integration. This method, 
supplemented with tree-level event simulations to compute the acceptances associated with the 
specific cuts of each experimental search, will allow us to explore efficiently the multi-dimensional 
parameter space of our model avoiding a time-consuming scan. Not all the production and decay 
processes that could be computed with this method are equally sizable, however. In sect. 3.2 we will 
present an estimate of the various processes based on the use of the Goldstone boson Equivalence 
Theorem [23], this will allow us to classify (in sect. 3.3) the channels which are more promising for 
the search of the top partners at the LHC. 

3.1 Production and Decay 

Given that the partners are colored they can be produced in pairs through the QCD interactions. 
The pair production cross-section is universal for all the partners and it can be parametrized by a 
function 

O-pairl^n-x) , (3.1) 

which depends uniquely on the partner's mass mx, for which we have analytical formulae. We have 
constructed Upair by interpolation using the HATHOR code [24] which incorporates perturbative 
QCD corrections up to NNLO. The values of the cross-section used in the fit are reported in Table 2 
for the LHC at 7 and 8 TeV center of mass energy. In this and all the other simulations we adopted 
the set of parton distribution functions MSTW2008 [25]. 

The other relevant process is the single production of the top partners in association with either 
a top or a bottom quark. This originates, as depicted in Figure 3, from a virtual EW boson 
V = {W^, Z} emitted from a quark line which interacts with a gluon producing the top partner 
and one third-family anti-quark. The possible relevance of single production was first pointed out 
in Ref. [26] . The relevant couplings have the form 

gxti^XnttR + gxtr^^LftL + gxh^XhthL , (3.2) 

where X = {T, B, X2/3, X5/3, T} denotes generically any of the top partners. At each vertex the EW 
boson V is understood to be the one of appropriate electric charge. Notice that there is no vertex 
with the hn because the latter state is completely decoupled in our model, we expect this coupling 
to be negligible even in more complete constructions. 

It is important to outline that the couplings gxt^ , gxt^ ^i^d gxbL '^^^ be computed analytically in 
our models. They arise from the interactions reported in Appendix B after performing the rotation to 
the physical basis of mass eigenstates. Since the rotation matrices can be expressed in a closed form 
the explicit formulae for the couplings are straightforwardly derived. The result is rather involved 
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a[fb] @ 


NNLO 




pair production 


M [GeV] 


^/s = 7 TeV 


= 8 TeV 


400 


(0.920) 1.41 xlO^ 


(1.50) 2.30 xlO^ 


500 


(218) 330 


(378) 570 


600 


(61.0) 92.3 


(113) 170 


700 


(19.1) 29.0 


(37.9) 56.9 


800 


(6.47) 9.88 


(13.8) 20.8 


900 


(2.30) 3.55 


(5.33) 8.07 


1000 


(0.849) 1.33 


(2.14) 3.27 


1100 


(0.319) 0.507 


(0.888) 1.37 


1200 


(0.122) 0.196 


(0.375) 0.585 


1300 


(4.62) 7.60 xlO-2 


(0.160) 0.253 



Table 2: Cross sections for the NNLO pair production of heavy fermions at -^s = 7, 8 TeV (the LO values 
are in brackets), with HATHOR [24]. 




Figure 3: The single-production diagrams. 

and for this reason it will not be reported here, however it is easily implemented in a Mathematica 
package. 

The single production cross-sections are quadratic polynomials in the couplings, with coefficients 
that encapsulate the effect of the QCD interactions, the integration over the phase-space and the 
convolution with the parton distribution functions. These coefficients depend uniquely on the mass 
of the partner and can be computed by Monte Carlo integration. Once the latter are known we obtain 
semi-analytical formulae for the cross-sections. The production in association with the b is simply 
proportional to Oxb^ '^^lile the one with i would be, a priori, the sum of three terms proportional 
to gxt^, 9xtn ^^"^ 9xtL ' 9xtji which account, respectively, for the effect of the left-handed coupling, 
of the right-handed one and of the interference among the two. However in the limit of massless 
top quark, mt <^ mx , the processes mediated by the left-handed and by the right-handed couplings 
become physically distinguishable because the anti-top produced in association with X will have 
opposite chirality in the two cases. Therefore in the limit — )• the interference term can be 
neglected. Moreover, the coefficients of the gxt^"^ and gxtR^ terms will be equal because the QCD 
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a[pb] 
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a[ph] 
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single production of tB + tB 


single production of bT + bT 


M [GeV] 


y/s = 7 TeV 


^/s = 8 TeV 


= 7 TeV 


= 8 TeV 


400 


(2.70) 3.10 


(4.32) 4.92 


(32.49) 43.47 


(47.83) 61.43 


500 


(1.49) 1.80 


(2.50) 2.97 


(15.85) 20.44 


(24.10) 33.10 


600 


(0.858) 1.06 


(1.49) 1.84 


(8.53) 12.89 


(13.55) 18.80 


700 


(0.511) 0.637 


(0.928) 1.15 


(4.60) 6.70 


(7.92) 11.34 


800 


(0.313) 0.399 


(0.590) 0.745 


(2.82) 4.01 


(4.58) 7.22 


900 


(0.194) 0.250 


(0.377) 0.497 


(1.60) 2.50 


(2.89) 4.48 


1000 


(0.121) 0.160 


(0.246) 0.325 


(0.956) 1.636 


(1.81) 2.83 


1100 


(0.075) 0.103 


(0.164) 0.215 


(0.604) 0.980 


(1.181) 1.72 


1200 


(0.048) 0.066 


(0.107) 0.146 


(0.377) 0.586 


(0.726) 1.23. 


1300 


(0.031) 0.043 


(0.072) 0.098 


(0.234) 0.386 


(0.463) 0.731 



Table 3: Cross sections for the NLO single production of B and T for a unit coupling, at ^/s = 7, 8 TeV (the 
LO values are in brackets), with MCFM [29]. 



interactions are invariant under parity. Thus the cross-sections will be very simply parametrized as 



o-smg(V6) = ioxbL^ (^Vb{rnx) , 



(3.3) 



in terms of few functions avti^x) and aybi'iT^x)- The charge-conjugate processes, in which either 
X t oi X b are produced, can be parametrized in terms of a similar set of coefficient functions. The 
only difference is the charge of the virtual V emitted from the light quark line. We thus have 



O-sing(Vt) = (gxtL) +{9Xtii) (JYttimx), 
asingiXb) = {gxtL^ crytf,{mx) , 



(3.4) 



where V"^ denotes the charge conjugate of the vector boson V. A similar way of computing cross 
sections of the W — b fusion type of single- production was carried out in Ref. [27] where they adapted 
the fitting functions of Ref. [28] to non-SM couplings. 

One might question the validity of the zero top mass approximation which allowed us to neglect 
the interference and parametrize the cross-section as in eq.s (3.3) and (3.4). We might indeed 
generically expect relatively large corrections, of the order of mt/mx- However the corrections are 
much smaller in our case, we have checked that they are around 1% in most of the parameter space 
of our models. The reason is that the interference is further reduced in our case because the left- and 
right-handed couplings are never comparable, one of the two always dominates over the other. This 
enhances the leading term, gxt^ °^ 9xtR^ ™ comparison with the interference gxt^ ' gxtjf Moreover 
this implies that eq.s (3.3) and (3.4) could be further simplified, in the sum it would be enough 
to retain the term which is dominant in each case. We will show in the following section that the 
dominant coupling is gxtR in the case of the fourplet (models M45 and M4i4) and gxt^ in the case 
of the singlet (models MI5 and MI14). 
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Figure 4: In red dashed: the cross sections of pair production. In green and blue the single production of 
the T(in association with a b) and of the X5/3(in association with a t), respectively in model MI5 and M45. 
The point chosen in the parameter space is ^ = 0.2, ci = 1 and y — I. The value of C2 is fixed, at each value 
of Mip, in order to reproduce the top quark mass. 

It total, all the single-production processes are parameterized in terms of 5 universal coefficient 
functions (Tyy±t, crzt and cr-^^±f,. Notice that a possible azb vanishes because flavor-changing neutral 
couplings are forbidden in the charge —1/3 sector as explained in the previous section. As such, 
the single production of the B in association with a bottom quark does not take place. We have 
computed the coefficient functions a\Y±t and a\Y±b, including the QCD corrections up to NLO, 
using the MCFM code [29]. To illustrate the results, we report in Table 3 the single production 
cross-section with coupling set to unity, for different values of the heavy fermion mass, and for the 7 
and 8 TeV LHC. The values in the table correspond to the sum of the cross sections for producing 
the heavy fermion and its antiparticle, on the left side we show the results for t B production, on 
the right one we consider the case of hT. In our parametrization of eq.s (3.3) and (3.4) the cross- 
sections in the table correspond respectively to (Jy/+t + Ciy-t and to (Jy^+i, + cTw-b- We see that the 
production with the h is one order of magnitude larger than the one with the t, this is not surprising 
because the t production has a higher kinematical threshold and therefore it is suppressed by the 
steep fall of the partonic luminosities. The values in the table do not yet correspond to the physical 
single-production cross-sections, they must still be multiplied by the appropriate couplings. 

The last coefficient function azt cannot be computed in MCFM and therefore to extract it 
we used a LO cross section computed with MadGraph 5 [30] using the model files produced 
with FeynRules package [31]. To account for QCD corrections in this case we used the k- factors 
computed with MCFM for the t B production process. 

In order to quantify the importance of single production we plot in figure 4 the cross-sections for 
the various production mechanisms in our models as a function of the mass of the partners and for 
a typical choice of parameters. We see that the single production rate can be very sizeable and that 
it dominates over the QCD pair production already at moderately high mass. This is again due to 
the more favorable lower kinematical threshold, as carefully discussed in Ref. [16]. 

Let us finally discuss the decays of the top partners. The main channels are two-body decays 
to vector bosons and third-family quarks, mediated by the couplings in eq. (3.2). For the partners 
of charge 2/3 and —1/3 also the decay to the Higgs boson is allowed, and competitive with the 
others in some cases. This originates from the interactions of the partners with the Higgs reported 
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in Appendix B, after the rotation to the physical basis of mass eigenstates. The relevant couplings 
can be computed analytically similarly to the Qti^ rX and gti^x- Thus we easily obtain analytical 
tree-level expressions for the partial widths and eventually for the branching fractions. In principle 
cascade decays X — )• X'V or X'H are also allowed, however these are never sizable in our model as 
we will discuss in sect. 3.3. 



3.2 Couplings to Goldstone Bosons 

Let us now turn to classify the relative importance of the various production mechanisms and decay 
channels described in the previous section. Since the partners are much heavier than the EW bosons, 
nix ^ nT-Wj their dynamics is conveniently studied by using the Equivalence Theorem, which applies 
at energies E ^ mw. To this end, we will momentarily abandon the unitary gauge and describe 
our model in the i?^-gauge where the Goldstone degrees of freedom associated with the unphysical 
Higgs components are reintroduced. The Higgs field is now parameterized as ^ 




The Equivalence Theorem states that, at high energies, the longitudinal components of the 
and of the Z bosons are described, respectively, by the charged and the neutral Goldstone fields (j)^ 
and (j)^. The transverse polarizations are instead well described by vector fields Wj^ and Z^, in the 
absence of symmetry breaking. However the transverse components give a negligible contribution 
to our processes, and this is for two reasons. First, their interactions emerge from the SM covariant 
derivatives and therefore these are proportional to the EW couplings g or g' . We will see below that 
the couplings of the longitudinal, i.e. of the Goldstones, are typically larger than that. Second, the 
transverse components can not mediate, before EWSB, any transition between particles in different 
multiplets of the gauge group. Indeed the couplings of the and fields are completely fixed 
by gauge invariance and therefore they are diagonal in flavor space. Only after EWSB do states 
from different multiplets mix and flavor-changing couplings like in eq. (3.2) arise. Therefore these 
effects must be suppressed by a power of e = {h)/f. This means that the transverse gauge bosons 
basically do not participate to the production and decay of the top partners: the decay will mostly 
be to longitudinally polarized vectors, while the virtual V exchanged in single production diagram 
will be dominantly longitudinally polarized. 

For our purposes, we can thus simply ignore the vector flelds and concentrate on the Goldstones. 
In the models with the fourplet, M45 (2.7) and M4i4 (2.11), the flrst source of Goldstone couplings is 
the term i ci (^_r) ^ 4* tR- One would naively expect this interaction to be the dominant one because 
it originates entirely from the strong sector without paying any insertion of the elementary-composite 
coupling y. Before EWSB the couplings are 



I — - — 



+ h.c. . 

(3.6) 



Notice that the Goldstone fields (j) ' in eq. (3.5) are not canonically normalized. Indeed the non-linearities in the 
Higgs kinetic term of eq. (A. 21) lead to a kinetic coefficient equal to sine/e, with e = (h)/ f. However this is irrelevant 
for the purpose of the present discussion. 
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It is not difficult to check that the interactions above respect not only the SM but also the full S0(4) 
symmetry of the strong sector. Eq. (3.6) contains derivative operators, therefore it is not yet suited 
to read out the actual strength of the interactions. However it can be simplified, provided we work 
at the tree-level order, by making use of the equations of motion of the fermion fields. After 
integrating by parts and neglecting the top mass, we find 



V2, 



ci 



/ 



V2 



- 1 TtR + ruBf, 



BtR + mxr 



2/3 



1\ 



mx, 



-5/3 ' 



+ h.c. , 
(3.7) 

showing that the strength of the interaction is controlled by the masses of the heavy fermions. 
Neglecting the elementary-composite coupling y, the masses all equal Mq,, and the coupling, modulo 
an 0(1) coefficient, is given by = M^/ f, as anticipated in the previous section. Once again we 
remark that this feature follows from the Goldstone boson nature of the Higgs. Indeed if the Higgs 
were a generic resonance, not a Goldstone, then it could more plausibly have a Yukawa g^:'^^IlitR 
vertex with strength dictated by the strong sector coupling . 

Those of eq. (3.7) are the complete Goldstone interactions in the limit of a negligible elementary- 
composite coupling y. However we can not rely on this approximation because we will often be 
interested in relatively light top partners, with g^ii < y ~ yt- It is straightforward to incorporate the 
effect of y, due to the mixing terms in eq.s (2.7) and (2.11) for model M45 and M4i4, respectively. 
After diagonalizing the mass-matrix, again neglecting EWSB, the Goldstone interactions for both 
models become 





M45, M4i4 


<^^-'^5/3L*iJ 

{p - i(tP)TL tR 
^-Bl tR 


V^cig^ 
cig^i, 



(3.8) 



which reduces to eq. (3.7) for y ^ g^. Notice that eq. (3.8) only contains couplings with the 
right-handed top quark. This is not surprising because the top partners live in SM doublets and 
therefore their only allowed Yukawa- like interactions are with the tR singlet. The couplings with 
the qi doublet emerge only after EWSB and are suppressed by one power of e. Therefore they 
typically do not play a mayor role in the phenomenology. Obviously the SM symmetry is respected 
in eq. (3.8), this explains the \/2 suppression of the X2/3 and of the T couplings compared with the 
ones of the X5/3 and of the B. 

The situation is different in the models with the singlet, MI5 and MI14 (2.12). In that case there 
is no direct contribution from the strong sector to the Goldstone coupling and all the interactions 



^"When considering a perturbation described by a small parameter to a Lagrangian, the use of the equations of 
motion of the unperturbed theory is equivalent to permorming field redefinitions of the form <3? — > $ + rjF[<^, d\. For 
example, to deal with the first term of eq. (3.6), the relevant redefinition is 



tR 



V2c 
J 



tR — —^hdTR 



This eliminates the derivative interaction and makes the first term of eq. (3.7) appear, ft also leads to new interactions 
with more fields that however are irrelevant for our processes at the tree-level. 
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are mediated by y. The couplings are 





MI5, MI14 


{p + i(tP)TRtL 


y 

y 



The top partner T now is in a SM singlet, therefore the interactions allowed before EWSB are the 
ones with the left-handed doublet. The \/2 suppression of the coupling with the top is due, once 
again, to the SM symmetry. One important implication of eq. (3.9) is that the T, contrary to the 
partners in the fourplet, can be copiously produced singly in association with a bottom quark. We 
will discuss this and other features of our models in the following section. 



3.3 The Most Relevant Channels 

We discuss here the most relevant production and decay processes of each top partner, identifying 
the best channels where these particles should be looked for at the LHC. Obviously one would need 
an analysis of the backgrounds to design concrete experimental searches for these promising channels 
and to establish their practical observability. We leave this to future work and limit ourselves to 
study, in section 4, the constraints on the top partners that can be inferred from presently available 
LHC searches of similar particles 

Let us first consider the models M45, M4i4 and analyze separately each of the new fermions. 

X5/3, together with X2/3, is the lightest top partner, it is therefore the easiest to produce. 
Production can occur in pair, via QCD interactions, or in association with a top quark through 
its coupling with a top and a W~^. The coupling, see eq. (3.8), is controlled by = mx^/g//, 
which grows with mass at fixed /. We thus expect single production to play an important 
role at high mass, where it is enhanced with respect to pair production by both kinematics 
and a larger coupling (at fixed /). This is confirmed, for a particular but typical choice of 
parameters, by the plot in Figure 4. 

Since it is the lightest partner, X5/3 decays to W~^t with unit branching ratio. The relevant 
channel for its observation is X5/3 — )• tW in association with a second top quark of opposite 
charge. The latter is present in both single and pair production processes. This results in clean 
signals consisting of either same-sign dileptons or trileptons plus jets. In the following section 
we will recast the LHC searches for these signals and obtain a limit on production. In 
addition to two top quarks and a W, pair production also leads to a second hard W while single 
production (see Figure (3)) features a light-quark jet associated with virtual W emission. 

Notice that the light-quark jet in single production is typically forward with a pr ^ "t-vf 
because the emission of the virtual W is enhanced in this kinematical region [16] . In practice 
this jet has the same features of the "tag jets" in VBF Higgs production and in VFVl^-scattering. 
The events are thus characterized by a forward isolated jet in one of the hemispheres. The 
relevant kinematical distributions are shown in Figure (5) for the production of a 600 GeV 
partner. Like in VBF or VKVF-scattering, one might hope to employ the forward jet as a tag 
to discriminate single production form the background. Ref. [16] argued that the main source 
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Figure 5: pr — V and energy distributions of the forward jets produced in a single production of the top 
partner with a mass 600 GeV. 

of forward jets in the background, QCD initial state radiation, tends to produce more central 
and less energetic jets, however further investigations are needed. Present LHC searches are 
designed for pair- rather than for single-production. Because of the r]^^^ and cuts that they 
adopt, they are thus weakly sensitivity to forward jets. We believe that it would be worth to 
explore the possible relevance of forward jets in designing the searches for top partners. 

X2/3 is also light and therefore easier to produce than the heavier partners. At the leading order, 
as eq. (3.8) shows, it couples with strength cig^ to the Higgs and Z bosons. The dominant 
decay channels are thus X2/3 — )• Zt and X2/3 — )• ht and BR{X2/3 Z t) ^ BR(X2/3 ht) ^ 
0.5. In model M45 the coupling to Wb vanishes exactly, while in model M4i4 the coupling 
is non-zero but suppressed by e ~ v/f. The decay X2/3 — )• Wb is therefore typically sub- 
dominant and can become relevant only in a corner of parameter space characterized by low 
mass, ye = 0(1) and ci < 1. Given that X2/3 — )• ht is probably difficult to detect, the search 
for X2/3 must rely on the decay mode X2/3 — )■ Zt, with Z further decaying to charged leptons. 
An extra suppression from the small branching ratio must then be payed. This disfavors the 
X2/3 signal compared to that of X5/3 , for which the branching ratio needed to reach the leptonic 
final state is close to one. 

X2/3 is produced in pairs via QCD interactions and singly via the ZX2/3t coupling,. In the 
latter case a top quark is produced in association. Both production modes lead to a resonant 
X2/3 — )• Zt plus one top of opposite charge. In the case of single production there will be a 
forward jet, as previously discussed in the case of X5/3. In the case of pair production there will 
be either a Higgs or a Z from the other partner. Another possible single production mode, in 
association with a b quark rather than a t, is strictly forbidden in model M45 and is suppressed 
by the small coupling to Wb in model M4i4. However single production in association with 
a 6 is kinematically favored over that with t. Kinematics then compensates the suppressed 
coupling and makes the two rates typically comparable in model M4i4, as shown in Fig. 6. 
By comparing with Fig. 4, we see that, in the case of X2/3, single production in association 
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Figure 6: Cross sections of the X2/3 pair (red dashed hne) and single production in association with a t (blue 
line) and with a b (green line) for the parameters choice: y = 1, ci — 1, (, — 0.2 in the model M4i4. 

with a f is suppressed compared to the case of X^^/^^^. This is mainly due to the \/2 factor in 
charged current versus neutral current vertices, see eq. (3.8). Moreover, the difference between 
the W and Z couplings, taking into account u- and the d-type valence quark content of the 
proton, further enhances by a ~ 1.2 factor the virtual W emission rate with respect to the Z 
rate. Combining this enhancement with the factor of 2 in the squared coupling, one explains 
the relative sizes of the X5/3 and X2/3 production cross sections. 

• T 

T is systematically heavier than X2/3, but the phenomenology is very similar. Therefore it 
will merely give a subdominant contribution to the X2/3 channels described in the previous 
paragraph. Indeed, by eq. (3.8), also T couples at leading order with equal strength to the 
Higgs and to the Z, leading to BR(r Zt) ^ BR(T ht) ^ 0.5. The coupling to Wb 
arises at order e, and it can be relevant, as explained for X2/3 above, thanks to the favorable 
kinematics of associated production with a b. 

One may in principle consider chain decays seeded by T — )• X2/3Z, T — )• X2/3/1 or T — X5/3W, 
given these channels are normally kinematically open. However the corresponding couplings 
are generically smaller than those controlling the direct decays to i/j. This is a straightforward 
consequence of the equivalence theorem and of SU (2) selection rules. The decays to tn, involve 
longitudinally polarized vectors and h, living in the linear Higgs doublet H: given the top 
partners are SU{2) doublets and tR is a singlet, the coupling respects SU{2) and so it arises 
at zeroth order in e. On the other hand, the transitions among top partners living in different 
SU (2) doublets obviously require an extra insertion of the Higgs vacuum expectation value. 
The resulting amplitudes are therefore suppressed by one power of e and the corresponding 
branching ratios negligible. 

• B 

B is even heavier than T, though the mass difference, niB — niT ~ y'^v^ /Atjib (see eq. (2.20)), 
is typically rather small. The most relevant decay mode is i? — )• Wt, mediated by the coupling 

^^Even though the two plots correspond to different models the couplings of the and X2/3 do not differ at 

leading order in models M45 and M4i4. 
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~ cigq, in eq. (3.8). Like in the case of T, SU{2) selection rules suppress the decay to WX2/s- 
Moreover, the decay B — )• WT, when kinematically allowed, proceeds either via a transverse 
W, with SM gauge coupling g < gq,, or via a longitudinal W, with effective coupling suppressed 
by e. Therefore also this decay is significantly suppressed. The decay i? — )• Z6 is forbidden 
because, as we explained in sect. 2.2.2, flavor-changing neutral couplings are absent in the 
charge —1/3 sector. The B ^ hb channel is forbidden in model M45 and suppressed by e in 
model M4i4. In the latter model it can play a role, but only in a corner of the parameter 
space. 

Single production, since the ZBb vertex is absent, is always accompanied by a top quark. The 
signature of single B production is therefore a resonant B — t- Wt plus an opposite charge top, 
the same final states of single X5/3 production. In the end, B production, single and pair, has 
the same signatures as X5/3 production: same sign leptons or trileptons plus jets. 

Let us now switch to models MI5 and MI14, where the only new heavy fermion is the T. 

• f 

T has a very rich phenomenology because it can be copiously produced through all the three 
mechanisms described above. We see in eq. (3.9) that T couples to both Zt and Wb, with a 
coupling of order y ^ yt/c2. It can therefore be singly produced either in association with a 
top or with a bottom quark. Notice that in the range C2 ~ 1 suggested by power counting, the 
trilinear coupling is of order yt, which is expected to be generically smaller than the strong 
sector coupling g^ that controls the single production of top partners in a (2,2). The bands 
in the left panel of Fig. 7, indicate the single prooduction cross section^^ for 0.5 < C2 < 2: 
comparing the blue band to the corresponding case of X2/3t and X^/^t production in models 
M45 and M4i4 , one notices, as expected, a typically smaller rate for models MI5 and MI14. 
While y ~ yt (c2 ~ 1) is favored by naive power counting, one can entertain the possibility of 
choosing y > yt (c2 < 1), for which the single production rate can be sizeable. However, for 
a given value of rrij, and /, there is a mathematical upper bound ymax on y determined by 
eqs. (2.28). The right plot in Fig. 7 shows that ymax grows with nij, and that it is comparable 
in model MI5 and model Ml 14. In the left panel of Fig. 7, the green line and the blue line 
shows, respectively for Tb and Tt, the maximal allowed cross section, which basically coincides 
with the choice y = ymax For such maximal values the single production cross section can 
be quite sizeable. 

Single production of a T-like partner was considered in the context of Little Higgs models 
in Refs. [32, 33], and more recently for composite Higgs models in Ref. [34], where it was 
also considered the possibility of using a forward jet tag as a handle for this kind of searches. 
The total cross section in this channel is favored over single production with a t by both 
kinematics and by the V2 factor in charged current transitions. Indeed, as shown in Fig. (7) 
associated Tb production dominates even over pair production in all the relevant mass-range 

^■^By fixing mt, £,, C2 and nij, the result for model MI14 and MI5 coincide. Indeed, by comparing the la- 
grangians (B.6) and (B.7), one notices that the gauge vertices and the mass spectrum of model MI14 equal those 
of model MI5 when the equality y^^'' sine = ^^114 gjjj2e/2 holds. 

^■^Note that, for a given m^, ymax does indeed correspond to the maximal value of the WfeT-coupling, while the 
coincidence is not exact in the case of the ZiT-coupling. 



24 



[GeV] 



[GeV] 



Figure 7: Left panel: cross sections for tlie different production mechanisms of T for the models MI5 and 
MI14 for ^ — 0.2. Red dashed: pair production; green line: Tb production with the maximal allowed coupling, 
green band: Tb production for 0.5 < C2 < 2; blue line: Tt production for the maximal allowed coupling, blue 
band: Tt production for 0.5 < C2 < 2. Right panel: maximal allowed y for the models MI5 (in yellow) and 
MI14 (in red). 



while single production with the t is rather smaU. The role of kinematics is especially important 
in this result, as the large Tb cross section is dominated by the emission of a soft 6, with 
energy in the tens of GeV, a regime obviously unattainable in the similar process wih a t. 
Indeed by performing a hard cut of order rrit on the pT of the h, the Th cross section would 
become comparable to that for Tt. Unfortunately the current LHC searches do not exploit the 
large inclusive rate of production with the h quark because they are designed to detect pair 
production. We will show in the following section that the acceptance of single production, 
with the cuts presently adopted is extremely low. We believe there is space for substantial 
improvement in the search strategy. 

Also concerning decays, all the possible channels are important in the case of T. It decays to 
Wh, Zt and ht at zeroth order in e, with a fixed ratio of couplings. By looking at eq. (3.9) 
we obtain BR(r ^ Zt) ^ BR(r ht) ^ 5 BR(r W b) ^ 0.25. Actually the branching 
fraction to Wb is even further enhanced by the larger phase space, though this is only relevant 
for low values of m^. Given that the branching fraction is larger, ideally the resonant Wb 
production would be the best channel to detect the T. However one should manage to design 
a search strategy to reject the background while retaining the signal. In particular one should 
retain as much as possible the contribution from the large single production in association with 
the b. A possibly cleaner decay channel could then be T — t- Zt with leptonic Z. 



4 LHC Bounds 

In this section we derive bounds on our models using the presently available LHC searches. Given 
that the top partners are heavy fermions coupled to top and bottom, we focus on the experimental 
searches for 4*'^ family quarks, which present a somewhat similar phenomenology We will make 

^^Significant bounds on the top partners could also emerge from unrelated studies like the searches of SUSY performed 
with the "razor" variable [35]. We thank M. Pierini for suggesting this possibility, obviously this is an interesting 
direction to explore. 
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use of the following searches for 4*^^ family quarks performed by CMS: 1) b' — )■ Wt with same-sign 
dileptons or trileptons in the final state [36]; 2) t' — )• Zt with trileptons in the final state [37]; 3) 
t' — )• Wb with two leptons in the final state [38] . 

In what follows, we quantify the impact of these three searches on our models, by adopting 
the following strategy. We compute separately the production cross-sections of the top partners, 
the branching fractions into the relevant channels and the efficiencies associated with the selection 
cuts performed in each experimental search. The cross-sections and the branching fractions at each 
point of the parameter space are encapsulated in semi-analytical formulae as described in section 3. 
The efficiencies are instead obtained numerically through parton level Monte Carlo simulations. 
Not having at our disposal a reliable tool to estimate the response of the detector, a fully realistic 
simulation of the hadronic final states would not be useful. Therefore we decided not to include 
showering and hadronization effects in our analysis. We partially took these effects into account 
by including the identification/reconstruction and trigger efficiencies reported in the experimental 
papers. 

4.1 Search for b' ^ Wt 

This search applies only to models M45, M4i4. The analysis of Ref. [36] aims at studying a 4*^^ 
family b' that is pair produced by QCD and is assumed to decay to Wt with unit branching fraction. 
The search is performed in the final state with at least one tagged b-jet and either same-sign dileptons 
or trileptons (e or /i). Three or two additional jets are required, respectively, in the dilepton and 
trilepton channels. Apart from the usual isolation, hardness and centrality cuts for the jets and the 
leptons, a hard cut is required on the scalar sum of the transverse momenta of the reconstructed 
object and missing pT- Ref. [36] reports the observed number of events in the two categories and 
the expected SM background. From these elements, given the efficiency of the signal in the two 
channels, one puts a bound on the pair production cross-section and eventually on the mass of the 
b'. With 4.9/6-1 of data at 7 TeV the bound is 611 GeV at 95% confidence level. Below we wih 
quantify the impact of this search on the parameter space of our models. 

The top partners contributing to the signal are X5/3 and B because, as shown in the previous 
section, they lead to two tops of opposite charge and to at least one extra W in both pair and 
single production. To derive the bound we must compute, for each partner and production mode, 
the efficiency of the signal in the dilepton and the trilepton channels as a function of the partner's 
mass. The total production cross-sections are computed semi-analytically at each point of the 
parameter space. Combining the cross-sections with the efficiencies we obtain the signal yield in 
the two channels that must be compared with the observed number of the events and with the 
expected background. We perform this comparison by computing the confidence level of exclusion 
(CL) defined through the CLs hypothesis test [39], as explained in some detail in Appendix C. At 
the practical level it is important that at the end of this procedure we obtain an analytical expression 
for the CL as a function of the fundamental parameters of our model. This makes very easy and 
fast to draw the exclusion bounds even if we work in a multi-dimensional parameter space. 

Efficiencies 

The first step is to simulate the signal processes. Rather than employing our complete model we 
have used a set of simplified MadGraph models containing the SM fields and interactions plus the 
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two relevant new particles - X5/3 and B - with the appropriate couplings to Wt responsible for 
the single production and the decay. We will employ the right-handed X^^Wtji or BWtpt vertices 
because, as we have shown in section 3.2, the top partners couple mainly to the tR. However, to 
make contact with Ref. [36], we simulated also the case of left-handed vertices because for a 4*^^ 
family b' the coupling originates from an off-diagonal entry of the generalized Vckm matrix and it is 
purely left-handed. We will see that the chirality of the couplings significantly affects the efficiencies. 

We generated parton level events without showering, hadronization and detector simulation. The 
events were analyzed using the cuts and the identification/reconstruction efficiencies for 6-tagging 
and leptons reported in [36]. We also included the trigger efficiency as an overall multiplicative 
factor. In order to account for the possible merging of soft or coUinear partons in a single jet we 
applied the anti-/cT clustering algorithm [40] for jet reconstruction. Not having enough information 
on how the r leptons were treated in the analysis we have accounted only for the missing energy from 
the tau decays while the jets and leptons candidates coming from taus were simply rejected. We 
checked that the inclusion of r-jets does not introduce appreciable differences, but r- leptons might 
affect our results. We have found that the most severe cut is the one on the transverse momenta of 
the leptons candidates, pT > 20 GeV. This is because most of the events which could contribute to 
di-(tri-) leptons contain exactly 2(3) charged lepton candidates, thus loosing only one of them causes 
the loss of the event. The number of generated jets per event is instead larger than the minimally 
required one and therefore the impact of the jet cut is less prominent. 

The signal efficiency is defined as the product of the cut efficiencies with the branching ratios of 
the t and the W to the required final states. The results are given in Tables 4 and 5 for different mass 
points. In these tables, the efficiencies of the pair-produced b' obtained in Ref. [36] are compared 
with the ones obtained for a left-handed coupling with our method. The accuracy of our simplified 
treatment of QCD radiation and detector effects is quantified by the level of agreement between 
these results. We see that the discrepancy is below 10% in the dilepton channel and around 30% 
in the case of trileptons. In view of these results we have decided to be conservative and to present 
our results by showing exclusion limits computed using our efficiency and also using an efficiency 
reduced respectively by 10% for dileptons and 30% for trileptons. From Tables 4 and 5 we also see 
that the efficiency in our model is significantly larger than the one for the 4*'* family b'. This is 
because the right-handed top (and the left-handed anti-top) produced in the decay in our models 
tends to produce more energetic charged leptons than a left-handed top. The lepton pT distribution 
is therefore harder and the cut pT > 20 GeV is more easily satisfied. Finally, we notice, somewhat 
surprisingly, that the efficiencies for the X5/3 and for the B partners are substantially identical. One 
would have expected some difference at least in the dilepton channel, since the two leptons come 
from the decay of a single heavy particle in the first case while they have a different origin in the 
second one. However this makes no difference in practice. 

Plots and results 

Now using the event analysis algorithm shortly described above we can compute the signal efficiencies 
for same-sign dileptons and trileptons in the framework of the model with a totally composite top 
right and top partners in a four-plet. For this we again employ simplified models with only two top 
partners but this time we use exact couplings corresponding to typical points in the parameter space. 
Therefore apart from the right-handed coupling which is still dominant there is a small admixture 
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M [GeV] 


X5/3 partner [%] 


B partner [%] 


4*1^ family b' [%] 


b' Ref. [36] [%] 


450 


1.90 ±0.05 


1.93 ±0.05 


1.65 ±0.04 


1.52 ±0.13 


550 


1.97 ±0.05 


1.98 ±0.05 


1.72 ±0.05 


1.71 ±0.14 


650 


1.96 ±0.05 


1.96 ±0.05 


1.85 ±0.05 


1.71 ±0.15 



Table 4: Efficiencies for the pair produced B and going to same-sign dileptons. Efficiencies contain 

the cuts losses, b-tagging performance and BR's of W boson. 



M [GeV] 


X5/3 partner [%] 


B partner [%] 


4*^ family b' [%] 


b' Ref. [36] [%] 


450 


0.88 ±0.02 


0.84 ±0.02 


0.69 ±0.02 


0.47 ±0.05 


550 


0.98 ±0.02 


0.94 ±0.02 


0.81 ±0.02 


0.56 ±0.05 


650 


1.04 ±0.03 


1.07 ±0.02 


0.82 ± 0.02 


0.63 ±0.06 



Table 5: Efficiencies for the pair produced B and X5/3 going to trileptons containing two opposite-sign 
leptons. Efficiencies contain the cuts losses, b-tagging performance and BR's of W boson. 

of the left-handed one. We present the results for X5/3 and B masses in a range 400 — 1000 GeV in 
tables 6 and 7 for pair and single production respectively. 

Now, by using the obtained efficiencies together with the method elaborated above for computing 
the cross sections, one can compute the number of signal events in dileptons and trileptons and check 
if it falls into the region allowed by Fig. 12. 

In Fig. 8 we show the excluded region in the {^,M-^5fi) plane, where = depending on 

whether the single production is suppressed (ci = 0.3) or enhanced (ci = 3) and whether also B 
contributes to the signal {Mb > -^^5/35 V — O-^) or not {Mb ^ Mxyg, y = 3). Fig. 9 shows the 
exclusion in terms of Mx^^^ and ci. Since, as was discussed in sect. 3.2, the leading contribution to 
single production couplings is the same for models M45 and M4i4, the excluded regions are also 
similar for both models. A difference shows up when ci ^ 1 and the hBb vertex of model M4i4 
becomes important thus decreasing BR(i? — t- Wt) and also when = 0{1) and higher order 
effects modify the single production couplings. The excluded regions are almost symmetric with 
respect to ci — )• — ci, which can be understood as follows. When only X5/3 production matters, the 
single production rate is proportional to |cip at lowest order in e. Higher order terms only matter 
in the region of small |ci| where the single production rate is anyway negligible and the bound is 
driven by pair production which is insensitive to ci. When B production matters, that is because 





dilept eff. 


trilept eff. 


M[GeV] 


for B [%] 


for X5/3 [%] 


for B [%] 


for X5/3 [%] 


400 


1.67 ±0.03 


1.61 ±0.04 


0.66 ±0.01 


0.67 ±0.02 


600 


1.96 ±0.03 


2.02 ±0.04 


0.93 ±0.01 


0.93 ±0.01 


800 


1.81 ±0.03 


1.86 ±0.04 


0.98 ±0.02 


0.97 ±0.02 


1000 


1.63 ±0.03 


1.63 ±0.04 


0.99 ±0.02 


0.96 ±0.02 



Table 6: Efficiencies for the pair produced B and going to dileptons and trileptons containing two 

opposite-sign leptons. Efficiencies contain the cuts losses, b-tagging performance and BR's of W boson. 
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dilept eff. 


trilept eff. 


M[GeV] 


for B [%] 


for X5/3 [%] 


for B [%] 


for X5/3 [%] 


400 


0.50 ±0.01 


0.49 ±0.01 


0.12 ±0.01 


0.12 ±0.01 


600 


0.68 ±0.01 


0.69 ±0.01 


0.22 ±0.01 


0.22 ±0.01 


800 


0.65 ±0.01 


0.74 ±0.01 


0.26 ±0.01 


0.25 ±0.01 


1000 


0.63 ±0.01 


0.70 ±0.01 


0.28 ±0.01 


0.27 ±0.01 



Table 7: Efficiencies for the single produced B and going to dileptons and trileptons containing two 

opposite-sign leptons. Efficiencies contain the cuts losses, b-tagging performance and BR's of W boson. 

niB — ^ "^^5/3' corresponding to y ^ g^. From eq. (3.8) it is then evident that in this regime 

the couplings of both particles are approximately oc ci , so that the signal yield is again symmetric 
under ci — ?■ — ci. The gray regions on the plots correspond to an estimate of the mentioned above 
error in the determination of the efficiencies. 

4.2 Search for t' -> Z t 

The search in Ref. [37] is designed to detect an up-type 4*^ generation quark t' pair-produced by 
QCD and decaying to Zt. The search is performed in the trilepton channel, with two same-flavor 
and opposite-charge leptons with an invariant mass around the Z pole. Moreover, at least two jets 
are required. Apart from the usual hardness and isolation cuts for jets and leptons an important 
event selection is performed with the variable Rt, defined as the scalar sum of the reconstructed 
momenta without including the two hardest leptons and the two hardest jets. Rt is required to be 
above 80 GeV. With 1.14 /^-^ of 7 TeV data the bound on the t' is of 475 GeV. 

All the top partners of charge 2/3 can contribute to this final state these are X2/3 and T 
in models M45 and M4i4 and T in models MI5 and MI14. Remember however that the masses 
and the couplings of X2/3 and of T are closely tight to those of, respectively, X5/3 and B. Namely, 
the masses are similar (or equal) and the couplings at the leading order (see eq. (3.8)) differ by 
a factor of \/2- Therefore the search of charge 2/3 states will constrain the same combinations 
of the fundamental parameters of the model. But the bound on the charge 2/3 partners which 
can be obtained using Ref. [37] is by far less stringent than the one from the b' — t- Wt search [36] 
described above. Given approximately 5 times less of analyzed data the limit on the production 
cross section of Ref. [37] is significantly looser than the one of the Ref. [36] . Moreover in our model 
the production yield of the charge-2/3 states is typically lower than the one of the A5/3 and of the 
B. This is because of the branching ratio suppression to reach the tZ final state and because the 
single production rate is smaller (see section 3.3). Therefore we can safely ignore Ref. [37] when 
constraining models M45 and M4i4. And moreover, for the reasons described above, we expect 
that, even by updating the search of Ref. [37] to the same integrated luminosity of Ref. [36], it would 
not become more important. 

Hence in the following we will only use the search for t' ^ Z t to constrain the parameters of 
models MI5 and MI14 . As in the previous section, we will obtain semi-analytical formulae for the 

Actually because of the quite loose cuts on the invariant mass of leptonic Z used in this search also the X^^^ and 
the B could contribute, however this effect is subdominant. 
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Figure 8: Excluded (95%CL) regions in the plane for the models M45 and M4i4, using the search 

for b' ^ Wt. In red: a = 0.3 and y = 3 {Mb > A^Xs/a), in blue: ci = 3 and y = 3 (Mb > Mx^,^), in 
green: ci = 3 and y = 0.3 (Mb ?J ^^5/3 for ^ > 0.1, Mb ^ for C ^ 0.1). Gray regions correspond 

to a variation of the dileptons and trileptons signal of approximately 10% and 30% respectively (see text for 
details). 

signal yield by computing the efficiencies at each mass point and multiplying with the production 
cross section computed in Section 3.1. Differently from the previous case, the search is performed 
in a single channel. Therefore we will not need any statistical analysis, we will just compare the 
computed signal yield with the 95% CL limit obtained in Ref. [37] which corresponds to 9.6 signal 
events. 

Efficiencies 

The efficiencies are computed with a MadGraph model which incorporates T and its couplings to 
Zt, ht and Wb. These are responsible for the two single production modes and for the decay. The 
T couples only to left-handed quarks, therefore in this case we will employ left-handed couplings to 
compute the efficiencies. Our results can thus be directly compared with the efficiencies reported in 
Ref. [37] because the coupling is left-handed also in the case of a 4^^ family quark. 

The T can contribute to the signal both in the pair and in the single production mode, provided 
that at least one T decays to Zt paying a branching fraction of around 1/4 (see Section 3.3). In the 
case of pair production all the decay modes of the second produced T {Zt, Wb or ht) are potentially 
relevant. We have computed separately the efficiencies in all three cases. The methodology of the 
analysis, and in particular the treatment of the experimental efficiencies, closely follows the one of 
the previous section. The results are shown in Table 8, our efficiencies contain the cut losses and the 
W, Z, h branching fractions to the required final state. The efficiencies listed in the ffist column of 
the table can be directly compared with the ones of Ref. [37], we have checked that the discrepancy 
is around 25% which corresponds to approximately 1.5a of the signal uncertainty obtained in the 
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Mx,„ [GeV] M^^^^ [GeV] 

Figure 9: Excluded (95%CL) regions in the {Mx^/^^ci) plane for ^ = 0.2 for the models M45 and M4i4, 
using the search for b' ^ Wt. In blue: y = 3 {Mg ^ Mx^^^), in green: y = 0.3 (M^ > Mxg/3). Black dashed 
lines correspond to the exclusions with ^ — 0.4. Gray regions correspond to a variation of the dileptons and 
trileptons signal of approximately 10% and 30% respectively (see text for details). 

Ref. [37]. 

We see in Table 8 that the efficiency for the single production with the b is extremely low, below 
1 %o. This is because the single production signal (see Figure 3) is characterized by three leptons 
plus one hard (6) jet from the top decay, plus one forward jet from the virtual W emission and a b 
from the gluon splitting. But the gluon splitting is enhanced in the collinear region, therefore the 
6-jet emitted from the gluon is also preferentially forward and with low px- In order for the event 
to pass the selection cut, that requires at least two jets with pT > 25 GeV and \rj\ < 2.4, at least 
one of the two preferentially forward jets must be central and hard enough, implying a significant 
reduction of the cross-section. However this is not yet the dominant effect, the main reduction of the 
signal is due to the cut Rt > 80 GeV discussed before. Indeed Rt is computed without including 
the two hardest leptons and the two hardest jets, which in our case means, since we have only 3 
leptons and typically only 2 jets, that the momentum of the softest lepton must be above 80 GeV. 
Therefore in the end the signal is completely killed. The situation is better for the single production 
with the t since one typically has more particles produced in this case and therefore the efficiencies 
are comparable with the ones of pair production. 

The situation is better for the single production with the t, the efficiencies are comparable with 
the ones of pair production (see Table 8). However, we have seen in section 3.3 (see fig. 7) that 
the rate of pair production is typically larger than the one of single production with the top, in the 
relevant mass range. Since the efficiencies are comparable we do not expect a sizable contribution 
from this process. The signal is totally dominated by the pair production and the BR(r Zt) is 
fixed to be about 1/4, as discussed in section 3.3. Therefore the bounds one can infer are mainly 
on nif, but a mild dependence on the other parameters (^ and y) is still residual in the BR. The 
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pair prod. eff. [%] 


single prod. efF. [%] 


M [GeV] 


Tf Zt Zi 


Tf Zt Wh 


Tf Zt hi 


Ttj 


Tbj 


300 


1.78 


1.22 


1.51 


1.13 


0.03 


350 


1.93 


1.47 


1.64 


1.17 


0.03 


450 


2.21 


1.81 


1.81 


1.25 


0.05 


550 


2.34 


1.93 


1.95 


1.30 


0.06 


650 


2.40 


2.12 


1.96 


1.35 


0.08 



Table 8: Cuts efficiencies for the charge 2/3 top partners going to trileptons for the case of pair production and 
different decay channels, and a single production for the cases of Z-t fusion(4*'' column) and W-b fusion(5*'' 
column). Efficiencies contain the cuts losses and BR's of W, Z and the SM Higgs boson. 

resulting bound is about nif > 320 GeV in both models MI5, MI14, and it is maximized at large ^ 
and small y rrij, > 350 GeV in model Ml 14. These bounds are not competitive with those coming 
from the t' ^ W b search, as we are going to discuss next. 

4.3 Search for f ^ 6 

The last experimental study that we are going to consider is the search for a 4*^* generation t' quark 
decaying to Wb [38]. The search is performed in the channel of two opposite sign leptons (away from 
the Z pole) with two tagged bottom quarks. A very important selection cut, which is needed to 
suppress the background from the top quark production, is that the invariant mass of all the lepton 
and 6-jet pairs, M/^, is above 170 GeV. This forbids that the lepton and the b originate from the 
decay of a top quark. Using data from 5fb^^ of integrated luminosity, a lower bound of 557 GeV 
was set on the t' mass [38]. 

In our models, only T can decay to Wb with a sizable branching fraction. We will therefore 
use Ref. [38] to put constraints on models MI5 and Mli4^®. The single production mode with 
the b is definitely not relevant in this case because it only leads to one lepton. The one with the 
t is also irrelevant because the second lepton would come from the decay of the top quark and it 
would not satisfy the cut on Mih. We are therefore left with pair production. Moreover, because of 
the Mfej > 170 GeV cut, and as was explicitly checked in Ref. [41], pair production contributes to 
the signal only if both T's decay to Wb. We are then left with the same channel, TT — )• WbWb, 
considered in Ref. [38]. The chirality of the coupling responsible for the decay is also the same as 
in the 4*'^ family case. Therefore the efficiencies can be extracted directly from Ref. [38] without 
any need for additional simulations. Given the efficiency and taking into account that the branching 
fraction of the T to Wb, we can easily compute the signal yield and compare it with the bound 
obtained in [38]. 

^^Also pair produced B and decaying to Wt contributes to the final states considered in Ref. [38]. However 

tlie resulting bound on these states is lower than the one obtained using Ref. [36]. In addition the signature used in 
Ref. [38] is insensitive to single production. Thus we do not expect any improvement of the bounds on the models 
M45 and M4i4 from this search. 
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Figure 10: Excluded (95%CL) regions in the (M^,^) plane, using the search t' -> Wb, for the models MI5 
and MI14 for y — 0.5 (green), y = 2 (blue) (corresponding approximately to C2 — 2 (green), C2 — 0.5 (blue)). 
In the gray dashed region there are no solutions for Mj,{y,^) when y = 2. 

Plots and results 

We show the excluded regions of the parameter space in terms of ^ and on the Fig. 10. The 
exclusion is stronger for larger y (and smaller C2) due to a larger BR[T — )• Wh) in this case. As was 
already discussed in the Section 3.2 the gauge interactions of the model Ml 14 are similar to the ones 
of the model MI5 and therefore the excluded regions are also similar. The difference is sizable in the 
region close to ^ = 0.5 where in the model MI14 interactions with a Higgs boson vanish according 
to Eq. (B.7) and therefore the BR of the competitive decay to Wh increases. The regions without 
solutions for T{y,^) when y is large correspond to those defined by the Eq. (2.28). 

Due to a larger amount of data analyzed and a higher BR of the T ^ W h decay mode the 
search of Ref. [38] gives a better constraint on the parameters of our models than the previously 
considered search T — )• Zt [37]. However one may expect that with increased amount of analyzed 
data the search for T — )• Z t can become competitive due to its sensitivity to single production. 

4.4 Summary of exclusions 

The results of the searches described above can be conveniently summarized by scanning over the 
values of the model parameters and selecting the most and the least stringent bounds on the top- 
partners' masses. The highest excluded masses of and correspond to the lowest value of y 
and highest ci and ^, and the opposite for the lowest exclusion. For T and B the highest exclusion 
corresponds to the highest y, ci and ^ and the opposite for the lowest exclusion. Maximal T mass 
exclusion is reached when y and ^ are maximal and the minimal exclusion is obtained for minimal y 
and ^. In Fig. 11, we show our results for the maximal and minimal exclusions obtained by varying 
the parameters in the ranges: y G [0.3,3], ci G [0.3,3] and ^ G [0.1,0.3]. 
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Figure 11: Maxmal and minimal bounds on the masses of top partners for y e [0.3,3], ci € [0.3,3] and 
^ G [0.1,0.3] for the models M45, MI5 (left pannel) and M4i4, MI14 (right pannel). Blue and green bars 
correspond respectively to high and low values of y. Black dashed lines correspond to the exclusions for the 
reference values ^ = 0.1, ci = 1, y = 1. 



5 Conclusions 

In this paper we described an approach to systematically construct the low-energy effective la- 
grangian for the lighest colored fermion multiplet related to the UV completion of the top quark 
sector: the top partner. Our construction is based on robust assumptions, as concerns symmetries, 
and on plausible assumptions, as concerns the dynamics. Our basic dynamical assumption, follow- 
ing Ref. [4], is that the electroweak symmetry breaking sector, or at least the fermionic sector, is 
broadly decribed by a coupling g(* and a mass scale m*. This assumption implies a well definite 
power counting rule. In particular the derivative expansion is controlled by inverse powers of m^,. 
In the technical limit where the top partner multiplet ^, is parametrically much lighter than the 
rest of the spectrum (m<i, ^ m^,), our power counting provides a weakly coupled effective lagrangian 
description of the phenomenology of ^ . The basic idea is that, in this case, the effects of the bulk of 
the unknown spectrum at the scale m=K can be systematically described by an expansion in powers of 
miii/m=K. The lagrangian obtained in this limit defines our simplified description of the top parters. 
One should however keep in mind that the most likely physical situation is one where —m\i, ~ m^j, 
where an effective lagrangian is formally inappropriate. In practice, however, we expect it to be more 
than adequate for a first semi-quantitative description of the phenomenology and certainly to assess 
experimental constraints. The comparison with explicit constructions supports this expectation. 

As concerns the symmetries of the strong sector, we considered the minimal composite Higgs 
based on the SO{5)/ S0{4:) coset. Furthermore we focussed on the simplest possibility where the 
right-handed top quark tR is itself a composite fermion. The leading source of breaking of SO (5) is 
thus identified with top quark Yukawa coupling yt. In our construction, we have fully exploited the 
selection rules obtained by treating yt as a small spurion with definite transformation properties. For 
instance the structure of the mass spectrum and the couplings are greatly constrained by symmetry 
and selection rules. In particular the pNGB nature of the Higgs doublet implies the couplings 
originating from the strong sector are purely derivative: at high energy, or for heavy on-shell fermions, 
these couplings are effectively quite sizeable and yet they do not affect the spectrum even accounting 
for (H) 7^ 0. If the Higgs were not treated as a pNGB a large trilinear would be associated with a 
large Yukawa coupling and the spectrum would necessarily be affected when (H) ^ 0. 
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Depending on the quantum numbers of the top partner multiplet ^ and of the composite operator 
O that seeds the top Yukawa in the microscopic theory, one can then consider a variety of models. 
We focussed on the four possibiUties shown in Table 1, which could be considered the simplest ones. 
Our method can however be directly applied to perhaps more exotic possibilities. For instance one 
exotic, but not implausible, case would he O = 142/3 with ^ in the symmetric traceless tensor of 
SO {4), that is ^' = 92/3. This case involves a top partner with electric charge 8/3, performing a 
spectacular chain decay to 3W~^ + b. 

Our effective lagrangian depends on a manageable number of parameters. Once the top mass is 
fixed, beside the Goldstone decay constant / and the partner mass m^, there remain, depending on 
the model, only one or two additional parameters. These parameters, ci^2, control the size of the 
trilinear couplings between ^, third family fermions, and vector bosons or Higgs. They thus control 
the decay and the single production of top partners. Moreover, naive power counting suggests a 
preferred 0(1) range for these parameters. This fact, coupled with the constraints due to symmetry, 
robustly implies a definite structure for the interactions vertices in each model. For instance, for 
the case where ^ spans an SO (A) quadruplet, the trilinear coupling to the Higgs doublet involves 
mostly a tji and is expected to be of the order of a strong sector coupling = m^/ f. Moreover it 
grows with making simple production even more important in the range of heavy ^. In the case 
of a singlet ^, the trilinear is of order yt and involves the left handed doublet (t, b)L- These details, 
including the chirality of the top, affect the collider phenomenology of the models and, consequently, 
the constraints from searches. 

Not only have we a few lagrangian parameters, but also they mainly affect phenomenology 
via their contribution to the trilinear couplings. Using this property we devised a semi-analytical 
way to efficiently simulate the contribution of single production to the signal. For any given mass 
miji, we numerically simulated single production once for all, assigning trilinear coupling equal 
to unity. The physical cross section was then obtained by folding this numerical result with the 
analytical dependence of the physical trilinear coupling on the model parameters. Thus, once the 
efficiencies associated with a given experimental search are known, the constraint in parameter space 
can be obtained analytically. We implemented the calculation of the cross-sections in a Mathematica 
notebook which is available on request. 

We applied our results to the presently available LHC searches. We focussed on the search for 4*^* 
family fermions, that have signatures similar to those of top partners, and recast them to constrain 
our models. The main results can be read from Fig. 9 and Fig. 10. The former figure shows that, in 
the relevant region ci = 0(1), single production has a mild but non-negligible impact on the bounds. 

In the course of our analysis, it became evident that there exists significant space for improvement 
in the search strategy if one wants to best constrain this class of models. The searches we used were 
tailored to pair production of heavy quarks, while single production of top partners has different 
features. First of all, in single production there is only one hard decaying object, with the t or b 
produced in association beeing often less hard (much more so in the case of a 6) : the cuts performed in 
present searches tends to penalize single production. Secondly, single production is always associated 
with a very forward jet, originating from the coUinear splitting of a typically valence quark into a 
longitudinally polarized vector boson. The resulting forward jet has exactly the same features of 
the tag jets of WW scattering. There is a good chance the use of the same tag in the searches 
for singly produced top partners would significantly extend the sensitivity. We also realized that 
the single production with the b is typically very large (see fig. 7) in the case of a singlet top 
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partner T, tagging this production mechanism would increase significantly the LHC sensitivity to 
this kind of particles. Ideally the best channel of detection would be the resonant Wb production 
from the T decay, accompanied by one forward jet from the longitudinal vector boson emission. 
The second b quark which is present in the reaction, which comes from the gluon splitting as in 
fig. 3, is typically quite soft both in p± and in energy. Thus it is probably strongly affected by 
QCD initial state radiation and difficult to detect. Building upon these considerations, it would be 
worth to undertake a thorough experimental analysis, including the effect of radiation and detector 
simulation, suitably designed for the search of singly produced top-partners. 

In the test of weak scale naturalness, the search for all possible fermionic top partners represents 
the other half of the sky. In this paper we have introduced a first systematic description of top 
partner phenomenology. The simplicity of the result should hopefully serve as a basis for future 
systematic experimental studies. As seen from our theorist's analysis, the present searches have 
already advanced well into the region suggested by naturalness. But there is no doubt somebody 
out there can do better. 
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A Explicit CCWZ construction for SO(5)/SO(4) 

Generators and Goldstone Matrix 

The generators of SO (5) in the fundamental representation are conveniently chosen to be 



L,R)IJ 



(A.l) 



(A.2) 



where ji (a 



1,2,3) are the S0(4) ~ SU(2)l x SV{2)r unbroken generators, while T' {i = 
1, . . . ,4) are the broken ones and parametrize the coset SO(5)/SO(4). An equivalent notation for 
unbroken generators which we will use is with a = 1, ... ,6. The indices IJ take the values 
1, . . . , 5. The normalization of the T'^'s is chosen as Tr[T'^, T^] = 6^^ . 

The and generators span respectively the SU(2)i and SU(2)/j subgroups, and obey the 
standard commutation relations 



pa rji/3 



(A.3) 
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The Tls are therefore identified as the generators of the SM S\J(2)l. Notice that in our parametriza- 
tion the unbroken T^'s are block-diagonal 



rpa 



e 




(A.4) 



and the generators obey the following commutation relation 

[r",r] = it'')jiT^ . (A.5) 

With these generators, the parametrization of the Goldstone boson matrix is explicitly given by 



U = U{U) = exp 



f ' 



nn^ 



L4x4 jJT" 



1 — COS ' 



n ''"^ / 



?) g™? 

cos J 



(A.6) 



where II'^ = 11*11. Under g € S0(5), the Goldstone matrix transforms as 

U{U) ^ U{U'^^^) = g-U{U)-h\U;g), 
where h{Il; g) is block-diagonal in our basis 



(A.7) 



(Ai 



with /i4 G S0(4). Under the unbroken S0(4) the II's transform linearly, using eq. (A.5) we get 
n* —7- (/14)* . Given our embedding of the SM group, the 11 four-plet can be rewritten as 



n 







n2 


1 




~ \/2 



(-i {hu - hl)\ 
i {hd - h\) 



(A.9) 



where 



H 



is the standard Higgs doublet of -1-1/2 Hypercharge. 
In the unitary gauge, in which 



hu = 0, 



hd 



h_ _ jh) + p 
~ V2 



(A.IO) 



(A.ll) 



where p is the canonically normalized physical Higgs field, the Goldstone boson matrix of eq. (A. 8) 
simplifies and becomes 

/I \ 

10 

[/= 1 . (A.12) 

cos^ sin^ 

VO -sin^ cos^/ 
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Given that we wiU have to gauge the SM subgroup of S0(5), we must consider also local trans- 
formations, g = g{x), in the above equation. We also have to define gauge sources 

= A^T^ ^ = g [A, + id^] g' , (A. 13) 

some of which we will eventually make dynamical while setting the others to zero. Explicitly, the 
dynamical part of A^ will be 

\ = -^W; [Tl + iTl) + ^W- [Tl - iTl) +g [c^Z^ + s^A^) T'l+g' {c^A^ - s^Z,) r| , (A.14) 

where Cw and Sw denote respectively the cosine and the sine of the weak mixing angle and g, g' are 
the SM couplings of SU(2)i and U(l)y. Notice that belongs to the unbroken S0(4) subalgebra, 
this will simplify the expression for the d and e symbols that we will give below. 

The d and e symbols 

Still treating A^ as a general element of the SO (5) algebra, we can define the d and e symbols as 
follows. Start from defining 

4 = 4^*) = U' [A^ + id^] U , (A. 15) 

this transforms under SO (5) in a peculiar way 

4 ^ Al^-^'-3'-3^ = 4^^) = h[A^ + id^] h' (A.16) 

Since h = h{Il; g) is an element of S0(4) as in eq. (A. 8), the shift term in the above equation, 
lives in the SO (4) subalgebra. Therefore, if we decompose A^ in broken and unbroken generators 

4 = -d;,r- e^lT", (A. 17) 

we have that transforms linearly (and in the fourplet of S0(4)) while the shift is entirely taken 
into account by e^. We have 

<^(/i4)Vm and e^ = ep'' ^ hi[ef,-id^,]h\. (A.18) 

Let us now restrict, for simplicity, to the case in which belongs to the S0(4) subalgebra, as 
for our dynamical fields in eq. (A.14). It is not difficult to write down an explicit formula for d and 
e, these are given by 

< =-A; + 4i "'°'^/^^' n't°v„n (A.19) 

where V^II is the "covariant derivative" of the IT field: 

v^rf = d^u' - iAi {t^y. w . (A.20) 

The first use we can make of the d^ symbol is to define the S0(5)-invariant kinetic Lagrangian 
for the Goldstone bosons, this is given by 

c^ = ^dX. (A.21) 
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In the unitary gauge of eq. (A. 11) and using eq. (A. 14) for the Goldstone Lagrangian becomes 



from which we can check that the field p is indeed canonically normalized and read the W and 

M 
/ 



Z masses mw = g/2f sin ^ , mz = mw/cw This fixes relation among {v) and the EW scale 



V = 246 GeV 

v = fsm^. (A.23) 

The symbol can instead be used to construct the CCWZ covariant derivatives, because the 
shift term in its transformation rule of eq. (A. 18) compensates for the shift of the ordinary derivative. 
Consider for instance the field ^I' defined in eq. (2.5) of the main text, which transforms in the 4 of 
S0(4), i.e. like ^' — )• /i4 • The covariant derivative is 

= df^^ + ie^^t"^. (A.24) 

The CP symmetry 

By looking at eq. (A) and remembering that CP acts as H{x) — )• H*{x^^^) on the Higgs doublet 
we immediately obtain the action of the CP transformation on the Goldstone fields 11 and on the 
Goldstone matrix U. It is 

n(x) ^ Q-nlx^^)), U{x) C75 • f/(x(^V a , (A.25) 

where C4 and C5 are respectively a 4 x 4 and a 5 x 5 diagonal matrices defined as 

C74 = diag(-1,+1,-1,+1), a = diag(-l,+l,-l,+l,+l). (A.26) 

In the above equations the superscript "(P)" denotes the action of ordinary spatial parity. Similarly, 
the ordinary action of CP on the SM gauge fields in eq. (A. 14) is recovered if we take 

^ C5- Ajf^ ■ C, . (A.27) 

Prom the above equations it is straightforward to derive the CP transformations of the d and e 
symbols defined in eq. (A. 17), 

cf^ ^ C^,{dPy , ^ C4 • (e(f )) • C4 . (A.28) 

In the fermionic sector, adopting for definiteness the Weyl basis, the CP transformation of the 
ql and of the tR are the usual ones 

X{x) ^ x^cn =i^o^2^*^^iP)^^ (A.29) 

for X = {^Lj^Lj^ij}- For the top partners, in the case in which they transform in the fourplet of 
S0(4) as in eq. (2.5), it is natural to define CP as 

^ C^i^fK (A.30) 

while for the case of the singlet we simply have ^ — t- vI'('^^). Notice that with this definition the 
charge eigenstate fields {T, i3, X2/3, X5/3} defined in eq. (2.5) have "ordinary" CP transformation as 
in eq. (A. 29); 
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B Fermion Couplings 

In this appendix we report the explicit form of the fermion coupHngs to gauge bosons and to the 
Higgs, in the unitary gauge defined by eq. (A. 11), that arise in our four models M45, M4i4, MI5 
and MI14, defined respectively in eq.s (2.7), (2.11) and (2.12). All the couplings are given before 
the rotations that diagonalize the mass- matrices. 

The first two terms, which are relevant for the models M45 and M4i4, are 



i "^n^itR = sin j{Xr^/r^)R\^^tR - sin jBrV^ Ir - sin jTr^Ir 

-IT— Sm'^{X2/3)R^tR + i [{X2/3)r - Tr] ^tR , 



(B.l) 



and the term with the symbol which we combine, for convenience, with the one from the covariant 
derivative in eq. (2.10) 



2 

3^^ 



.9 ~ 0^^)^5/3-^-^5/3 

C7;; Z O 



h 



+ — (^ cos ^ - '^sl)T$T + — (- J cos ^ - ^4)^2/3^^2/3 



+ 



9 

^/2 



cos 



h 



T + sin^ 



h 



X2/3 



+ "photon couplings" 



■ 2 ^ ^ , 2 ^ Y 



+ h.c. 



(B.2) 



The couplings to the photon are not reported explicitly in the above equation because they are 
simply the standard ones, being completely fixed by the U(l)em residual gauge symmetry. 

In addition to the ones in eq. (B.2), non-derivative couplings with the Higgs field emerge in model 
M45 from the terms 



yc2f(QlYUi5tR 



yf HBR + yftL 

2/C2/ h- 
— sm jtLtR . 



cos^ JjTr + sin^ ^(^2/3)7? 



(B.3) 



while in model M4i4 we have 



yf{Qiy'UiiUj5^'ji 

'-^(Qiyuj.Uj.tR 



yf cos fb,BR + ft, 
yc2f . 2h- 

-2Vf""7'"'^- 



(cos - + COS— )Tr 



(cOSy -COSy)(X2/3)R. 



(B.4) 



For the models with the singlet, MI5 and MI14, the only gauge-fermion interactions come from 
the covariant derivative 



3 Can 



slT^T + -eT4T . 



(B.5) 
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The couplings with the Higgs come instead from 

yc2f(QlYUi^tR -- 

for model Ml 5 and from 
for model M4i4. 

C Statistical tools 

In the analysis performed in the Ref. [36] the CLs method is used to obtain the exclusion confidence 
intervals for the mass of the b' quark. However this exclusion is made in terms of the pair production 
cross section assuming some fixed ratio between the yield in dileptons and trileptons channels. In our 
case this ratio depends on the relative strength of the single and pair production and can significantly 
deviate from the one used in the Ref. [36]. Thus we want to re-do part of the experimental analysis 
in order to extract a more model-independent exclusion in terms of the number of di- and trileptons 
separately. Though we are not restricted to using the CLg only, we think that this method is well 
suited for constraining the parameter space of our model. 

To use the CLg we first construct a test statistics g as a log-ratio of probability density for the 
signal+background hypothesis to the background hypothesis: 

i=2L3l ^ ' ^ 

where rii - number of observed in the pseudo-experiment di- and trilepton events, Si and bi - number 
of the signal and background events respectively. The distribution P for a small number of events can 
be taken as a Poissonian modified due to the presence of the uncertainties. The largest uncertainty 
in the experimental analysis comes from the background estimation and can be accounted for by 
taking a marginal probability density defined as 

P inM + = JvinM + ^A) InMiu., 1, 5.Jd.^ (C.2) 

where V stands for a Poissonian distribution of the observed number of events and InM for a 
log-normal distribution of the nuisance parameters centered at the value 1 with a variance corre- 
sponding to a relative error in the background estimation 5^, = S^Jbf. The analogous definition is 
taken for the background-only probability distribution P(ni\bi). 

The confidence level of the signal+background (backround only) hypothesis is defined as 

/•oo 

CL,t,(^t>) = / Psb(b){<l)dq (C.3) 

Jqobs 



yf . h ~ 
yc2f . h 



(B.6) 



yf . 2h ~ 
yc2f . 2h 



(B.7) 
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Figure 12: Excluded with 95%CL values of the signal in same-sign dileptons and trileptons channels (gray 
area) and the maximal allowed values of the di- and trileptons yields according to Ref. [36] (green points). 

where Psb{b)iQ) is a probability density of q which corresponds to distributed according to the 
signal+background (background only) hypothesis and q°^'^ corresponds to the observed number of 
events n°^^. Finally the exclusion confidence level for the signal Sj is: 

CL^xci = 1-5^ (C.4) 

Obtained in this way confidence intervals coincide with those given in the Ref. [36] with a relative 
deviation of excluded pair-production cross section less than 5%. The difference can be caused by 
our simplified treatment of the nuisance parameters, i.e. neglecting the signal uncertainties and 
assuming that the backgrounds of di- and trileptons are completely uncorrelated. 

Using the given above definition we find a region in a plane {32,83) excluded with 95%CL 
(Fig. 12). 
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